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m ■ ABSTRACT 

o : 

We have developed a new stellar population synthesis model designed to study early- 
type galaxies. It provides optical and near-infrared colors, and line indices for 25 absorp- 
tion lines. It can synthesize single age, single metallicity stellar populations or follow the 
galaxy through its evolution from an initial gas cloud to the present time. The model 
incorporates the new isochrones of the Padova group and the latest stellar spectral li- 
braries. We have applied our model to new data for a set of three early-type galaxies, to 
find out whether these can be fitted using single-age old metal-rich stellar populations, 
as is normal practice when one uses other stellar models of this kind. The model is 
extensively compared with previous ones in the literature to establish its accuracy as 
well as the accuracy of this kind of models in general. 

Using the evolutionary version of the model we find that we cannot fit the most 
metal-rich elliptical galaxies if we keep the IMF constant and do not allow infall of gas. 
We do however reproduce the results of Arimoto & Yoshii (1986) for the evolution of 
the gas, and produce colors, and, for the first time with this type of models, absorption 
line-strengths. It is in fact possible to fit the data for the elliptical galaxies by varying 
the IMF with time. Our numerical model is in good broad agreement with the analytical 
simple model. We prefer however to calculate the evolution of the gas numerically instead 
of using the simple model, since it offers more flexibility, and even improved insight, when 
comparing with observations. In the present paper we describe the model, and compare 
a few key observables with new data for three early-type standard galaxies. However 
the data, as well as our fits, will be discussed in much more detail in a second paper 
(Vazdekis et al. 1996), where some conclusions will be drawn about elliptical galaxies 
on the basis of this model. 



Subject headings: Elliptical Galaxies, Galactic Evolution, Stellar Evolution, Stellar Spec- 
troscopy, Chemical Evolution, Metallicity, Spectral Energy Distribution 
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1. INTRODUCTION 
1.1. Population Synthesis 

The information that is available about galaxies 
primarily relates to their morphology, internal kine- 
matics, and spectral energy distribution. The further 
away one goes, the less important become the first 
two, compared to the third. Clearly, if one wants to 
study galaxy evolution, study of the electromagnetic 
spectrum is of maximum importance. Since the spec- 
trum of a galaxy generally consists of a combination of 
stellar spectra, emission from gas, and possibly non- 
thermal radiation, partly extinguished by dust, and 
since also stellar spectra exist in a very large number 
of varieties, it is clear that understanding the spec- 
trum of a galaxy is very difficult. 

To study complicated spectra it is necessary to first 
understand the spectral energy distributions of rela- 
tively simple objects. For this reason we will discuss 
here a model that analyzes the spectra of early-type 
galaxies. These objects appear to contain relatively 
little dust extinction and gaseous interstellar medium, 
and to have little recent star formation. They have, 
for these reasons, been the most studied objects in 
the current literature on population synthesis. There 
have been a number of accepted ways to attack the 
problem of understanding the spectrum of an ellipti- 
cal. To understand why we adopt our current method, 
we will summarize shortly some of the population syn- 
thesis methods most often used in the literature. 

A stellar population synthesis program tries to find 
a combination of stars for which the integrated spec- 
trum agrees with the observed spectrum of the object 
under study. In practise the problem is often under- 
constrained, i.e. a number of combinations of stars 
can be found which are able to fit the spectrum. To 
overcome this problem one generally forces the solu- 
tion to obey certain constraints. These range from 
simple continuity requirements (e.g. the luminosity 
function should decrease monotonically) to the re- 
quirement that the distribution of stars is determined 
completely by stellar evolution calculations. Models 
with very few physical constraints are generally called 
empirical population synthesis models, as opposed to 
evolutionary models. 

Empirical models have been used with some suc- 
cess by Spinrad & Taylor (1971), Faber (1972), O'Con- 
ncll (1976, 1980) and Pickles (1985). These papers 
often make use of linear programming to obtain their 
results. Some workers, notably Bica (1988), have at- 
tempted to take into account evolutionary effects by 
using, as units of population, distributions of stars ob- 
served in clusters of our Galaxy, instead of individual 



stars. 

Evolutionary models use a theoretical isochrone 
or HR diagram, convert isochrone parameters to ob- 
served spectra in some way and, finally, integrate 
along the isochrone. They all need to make an as- 
sumption about which initial mass function IMF to 
use. Also, the models need a recipe prescribing when 
the stars have been formed. Since the IMF is not very 
well known at the present time, its treatment is not 
very different from one model to another. However, 
as far as the star formation rate (SFR) is concerned, 
some models assume that all stars are formed at the 
same time, others prescribe that the SFR has to de- 
crease exponentially with time, while still others ex- 
plicitly try to describe the whole formation of a galaxy 
from a gas cloud and form stars when the physical 
conditions in the gas are adequate. Examples of this 
kind of evolutionary models can be found in Tinsley 
(1968, 1972, 1978a,1978b,1980), Searle et al. (1973), 
Tinsley & Gunn (1976), Turnrose (1976), Whitford 
(1978), Larson & Tinsley (1978), Wu et al. (1980), 
VandenBerg (1983), Bruzual (1983,1992), Stetson & 
Harris (1988), Renzini & Buzzoni (1986), Rocca- 
Volmerange & Guiderdoni (1987,1988,1990), Guider- 
doni & Rocca-Volmcrange (1987,1990,1991), Yoshii & 
Takahara (1988), Rocca-Volmcrange (1989), Buzzoni 
(1989), Chariot & Bruzual (1991), Lacey et al. (1993) 
and Bruzual & Chariot (1993). Some of these models 
not only predict colors but also line-strengths notably 
those of Peletier (1989) and Worthey (1994). 

There are also models which combine evolutionary 
population predictions with considerations of chemi- 
cal evolution. These models follow the evolution of 
the gas and make use of isochrones of more than 
one metallicity (solar). Examples of these chemo- 
evolutionary population synthesis models are Arimoto 
& Yoshii (1986,1987), Casuso (1991), Bressan et al. 
(1994) as well as the model presented in this paper. 
Note that for this type of models only the global 
metallicity, Z, is normally taken into account to deter- 
mine the stellar populations. However it is in princi- 
ple possible to follow the abundance distribution sep- 
arately for each of several important elements in the 
gas. Examples of these so called chemical evolution 
models can be found in Larson (1972), Matteucci & 
Tornambe (1987), Tosi (1988) and there are many 
more. However, calculating colors and especially ab- 
sorption line-strengths in such models has not been 
attempted up to now. Finally, among this kind of 
models there are a few which combine chemical evo- 
lution with dynamics, e.g. Theis, Burkert & Hensler 
(1992). 

One cannot say that one type of model is better 
than another. In general, more observables and phys- 
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ical parameters can be calculated if one makes more 
assumptions. If no assumptions are made about the 
physics, as in the empirical models, one may end up 
with solutions that are unphysical. If however wrong 
assumptions arc made, one will not learn anything 
about the stellar evolution history either. We show 
later in this paper that our results can be reproduced 
using the simple analytical model and possibly with 
infall. This means that we could have replaced the 
part of the model that deals with the evolution of the 
gas by some analytical calculations. However, we have 
preferred to build up the numerical machinery, since 
it offers much more flexibility, and even improved in- 
sight. It is clear that our understanding of all these 
aspects is improving with time, which means that 
the models that are to be applied can legitemately 
be more and more complicated, and in this context 
we have developed the evolutionary population syn- 
thesis code presented in this paper. For the reasons 
mentioned above it should never be used as a static 
black box out from which a theoretical fit to the data 
is to be taken, but as an evolving tool, which might 
help in disentangling stellar populations in a compos- 
ite system. 

1.2. Evolutionary synthesis in general, and 
its problems 

If one wants to calculate the final spectrum of a 
galaxy that has evolved from a gas cloud, one has to 
integrate over time the spectra of all stars that are still 
living at the current time. The number of stars formed 
at a certain epoch is determined by the star forma- 
tion rate. Little is known about this, so many models 
give it a prescribed form, and let it decrease e.g. ex- 
ponentially In this paper we assume that the SFR is 
proportional to the gas density (Schmidt 1959). The 
gas density itself and the chemical evolution of the 
gas is calculated taking into account the original gas, 
and the metal-enriched gas that is ejected by stars. 
The yields for the various elements needed for this 
calculation are not especially well known, and bet- 
ter knowledge would significantly improve the model. 
Other factors that may affect the SFR, such as inflow 
or outflow, arc not considered in the context of the 
present models. 

The stars living at the current epoch contain stel- 
lar populations with a mixture of ages and metallic- 
ities. To calculate the final spectrum we decompose 
the stars into single stellar populations (SSP) each of 
a single age and mctallicity, and calculate their spec- 
tra. To do this, one needs in the first place theoretical 
isochrones. The parameters of the isochrones depend, 
amongst others, on opacities of ions and molecules, 
and are reasonably well known for the early stages 



of stellar evolution. However for later phases such 
as the AGB and the Post-AGB, evolutionary calcu- 
lations are very complicated, and could still benefit 
from significant improvement. The isochrones are 
much better for solar and sub-solar metallicity than 
for metal rich stars, because the former can be, and 
have been, tested observationally on globular clus- 
ters. In general, the relative composition of the el- 
ements heavier than Helium in the models is close to 
solar; only recently have people included for example 
oxygen-enhancement (Vandcnberg 1992) or a-elemcnt 
enhancement (Weiss et al. 1995). 

The following step is to obtain a spectrum for a 
star with physical parameters given by the isochrones. 
Model atmospheres needed for this (e.g. Kurucz 
1992) appear to be reasonably reliable in the blue. 
In the red, molecular opacities make them, for the 
time being, much less reliable. For that reason sev- 
eral authors (Faber et al. 1985, Gorgas et al. 1993 
and Worthey et al. 1994) have developed a method 
that depends much more on observations. They use 
a grid of observed stars with various theoretical pa- 
rameters to determine fitting functions that can be 
used to calculate an absorption line index for any 
combination of theoretical isochrone parameters. The 
problem however is that these fitting functions at the 
moment are available only for a limited number of ab- 
sorption lines, since one needs many stars of various 
types and metallicitics to make them. 

If one wants to carry out population synthesis, one 
needs to cover as large a wavelength range as possible, 
and also at as high a spectral resolution. This is to 
cover many colors and line indices, since in principle 
every color and absorption line is affected in a differ- 
ent way by the parameters above: the SFR, the IMF, 
the metallicity and abundance ratios etc. In practice 
a model with complete coverage from UV to near-IR 
is hard to implement, because of lack of fitting func- 
tions, color calibrations etc., and the observations are 
hard to obtain. In this paper we are concentrating on 
early-type galaxies, for which high spectral resolution 
is less important, due to their large velocity disper- 
sions. We have produced model output for colors and 
lines that are relatively easy to obtain and reproduce, 
and that allow us to separate effects due to various 
relevant physical parameters. To some extent we are 
limited by the availability of libraries of stars, espe- 
cially for the fitting functions, so this aspect too, can 
be significantly improved in the future. 

The layout of the paper is as follows. In Section 2 
we will describe how we obtain colors and absorption 
line indices for a single stellar population, emphasiz- 
ing differences from previous studies, and especially 
any improvements. In Section 3 we introduce our 
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chemical evolutionary model. In Section 4 we apply 
our model to a limited set of data from three standard 
galaxies, for which the data is presented in an accom- 
panying paper (Vazdekis et al. 1996, Paper II). We 
perform the fits for an SSP (the static option) and for 
a full evolutionary case. In Section 5 the conclusions 
are presented. 

Finally, in Paper II, apart from presenting the ob- 
servational basis of the data-set used here, we carry 
out a comprehensive fit of the whole set of indices and 
colors on the basis of the scenarios suggested in this 
paper. In particular, we compare the fits obtained 
assuming a single-age single-mctallicity stellar popu- 
lation or using the full chemical evolution model. We 
also discuss the relations between elements. 

2. A single-age, single-metallicity stellar pop- 
ulation 

The stellar content of a model galaxy calculated 
by our full evolutionary model consists of stars of 
various ages and metallicities. In practice its output 
spectrum is the integrated output of a large number 
of spectra of stellar populations each with a single 
age and metallicity (SSPs). In effect, the SSPs can 
be thought of as the building blocks of our stellar 
population model, although in an evolutionary model 
these SSPs are not independent but linked by physi- 
cal parameters governing star formation, which makes 
the evolutionary method intrinsically more meaning- 
ful. Many authors, e.g. Worthey (1994) (hereafter 
W94), restrict themselves to calculating a single SSP, 
which is less general, but whose simplicity can give 
insights in the stellar populations of a galaxy. We 
will refer to a single age single metallicity model as 
the static case of our stellar population model. In this 
section we will explain in detail how integrated colors 
and line strengths are calculated for this static case. 

This model works as follows. Assuming an age and 
a metallicity it calculates the present distribution of 
the stars according to the required isochrone (see next 
subsection) and weights the light of the stars on this 
isochrone with the assumed IMF (as defined in Sec- 
tion 2.2). Then it calculates the integrated colors and 
absorption lines, weighting each individual contribut- 
ing star by its luminosity (as explained in Sections 2.3 
and 2.4). We present the results of our model for a 
grid of ages and metallicities and compare them with 
the literature. 

2.1. The stellar data 



2.1.1. General 

Our models are based on stellar evolution theory 
and require a set of isochrones to predict the distribu- 
tion of the stars in the HR diagram at a given time. 
In this work we have adopted the large grids of the- 
oretical isochrones provided by Bertelli et al. (1994) 
(hereafter BBCFN). Their initial compositions have 
been chosen according to the empirical law of Pagel 
(1989) AY/AZ = 2.5 and arc [Z=0.0004Y=0.23], 
[Z=0.001,Y=0.23], [Z=0.004,Y=0.24], [Z=0.008,Y= 
0.25], [Z=0.02,Y=0.28] and [Z=0.05,Y=0.352]. All 
these sets of isochrones were computed using the most 
recent radiative opacities by Iglesias et al. (1992) 
(OPAL) except those of [Z = 0.001, Y = 0.23] which 
are based on models with the radiative opacities by 
Hucbncr et al. (1977) (LAOL). These isochrones 
cover a wide range of ages and of masses, from ap- 
proximately 4 x 10 6 to 17 x 10 9 yrs and 0.6 to - 72 M© 
respectively. Another main advantage in the applica- 
tion of these sets of isochrones is that they cover all 
the stellar evolutionary stages: from the ZAMS, be- 
yond the red giant tip until the white dwarf stage after 
the planetary nebula phase. For more details about 
these isochrones we refer the reader to BBCFN. 

The set of isochrones with Z = 0.001 can be safely 
used because the effect of varying from LAOL to 
OPAL opacities at low metallicities has been esti- 
mated to be small by Alongi et al. (1993) (see also 
BBCFN). 

2.1.2. Use of the isochrones of BBCFN 

In order to be used in our spectrophotometric stel- 
lar population synthesis models (see next sections), 
the isochrones must contain the following quantities: 
the age, the initial and the present mass along the 
isochrone, the effective temperature and the gravity. 

Because the electronic file with the isochrones pro- 
vided by the Padova group, instead of the initial mass 
contains the indefinite integral of the IMF by number 
over the mass (FLUM), as defined in BBCFN, we first 
obtain the initial mass by inverting the equation 

FLUM = — (1) 

1 — a 

where a — 2.35. Now, from the initial stellar mass 
and the mass at a given point on the isochrone, we 
obtain the gas mass fraction ejected by a star of mass 
to and metallicity Z until time t: R z (m,t) (see Sec- 
tion 3.2.1). 

Since the gravities are not in the electronic data 
base of the isochrones we have calculated them using: 

9=i^°G)^Tt ss (2) 
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where a and G are the Stcfan-Boltzmann and the 
gravitational constants respectively and mi SOC is the 
mass along the isochrone. 

2.1.3. The very low-mass stars 

As the low-mass limit of the isochrones used is 
0.6M©, we decided to complete them with stars of 
lower masses. The reason for this is that, even if 
these stars do not individually contribute substan- 
tially, their integrated light may not be negligible (see 
Table 5, above all in the near IR for old systems, in 
those cases where we require a population model with 
an IMF with high slope. 

To complete the lower-mass MS of the isochrones 
of composition Z=0.008, Z=0.02 and Z=0.05 we used 
the results for low-mass stars, for Z=0.004, Z=0.01, 
Z=0.02 and Z=0.03, generously provided to us by 
Pols (private communication). The corresponding 
tracks for solar composition are described in Pols et al. 
(1996). These tracks have been calculated using a re- 
newed version (Han, Podsiadlowski & Eggleton 1994) 
of the evolution program of Eggleton (1973), incorpo- 
rating an improved equation of state and making use 
of radiative opacity tables from Alexander & Ferguson 
(1994) for the low temperatures. Since the stars with 
masses below 0.6M Q are always on the ZAMS and 
show almost no evolution, it is straightforward to in- 
corporate them into the lower MS via the appropriate 
isochrones. However, for the isochrones with metallic- 
ity Z=0.008 we first linearly interpolated the low-mass 
stars between Z=0.004 and Z=0.01. Finally, to com- 
plete the low-MS of the super metal rich isochrones 
(Z=0.05) we were forced to extrapolate linearly from 
the low-mass stars of Z=0.02 and Z=0.03. As a test 
we see that in the theoretical HR-diagram there is 
good agreement between stars above (from BBCFN) 
and below 0.6M Q . 

2.2. The initial mass function. 

In this work we adopt two IMF shapes: 

i) Unimodal: 

$(m) = /3m~^ (3) 
where n for the solar neighborhood is equal to 1.35 
(VSaipeter), and (3 is a constant. 

ii) Bimodal: 

$(m) = (30A-» to < O.2M 

$(m) = (3p(m) 0.2 < to < O.6M (4) 

$(m) = /3m- 1* m > O.6M 

where p(m) is a spline for which we calculate the cor- 
responding coefficients solving the following system 
obtained by the boundary conditions: 




-1 o 1 

log m 



Fig. 1. — Plot of our unimodal and bimodal IMF with 
slope 1.35. Also represented are the data of Scalo 
(1986). 
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Fig. 2. — Plot of the bolometric correction BC(K) as a 
function of the color V — K for solar composition. The 
two dotted portions are linear interpolations to bridge 
discontinuities at V — K ~ 0.8 and V — K ~ 3, which 
result from using the different sources indicated. 
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p(0.2) = 0.4-** 
p'(0.2) = 
p(0.6) = 0.6-** 
p'(0.6) = -fi0.6- 



(5) 



This bimodal IMF is based on observational results 
of Scalo (1986) and Kroupa et al. (1993). 
Normalizing our IMF to unity we have 



/ 

J 771) 



l3m D t i dm + 



i 



/3m ^dm 



(6) 



where m; and m u are, respectively, the lower and up- 
per mass limits of the stellar range considered (in this 
paper we choose 0.0992M Q and 72 Mq for these lim- 
its). We obtain the constant (3 via numerical integra- 
tion. 

In Fig. 1 we plot our IMFs for \i = 1.35, comparing 
them with the data of Scalo (1986). 

2.3. The integrated colors 

The integrated flux of an SSP at a wavelength A, 
F\ TOT can be calculated using: 



Fx TC 



J m L 



N(m,T G )F x (m,T G )dm (7) 



where itit g is the mass of a star whose lifetime is Tq, 
F\{m, Tq) the flux of a star with mass m and age Tq, 
and N(m, Tq) is the number of these stars. This final 
distribution of stars is calculated via 



N(m,T G ) 



$(m) 



m 



M G Am 



(8) 



where M G is the mass of the galaxy (or of a given 
region of the galaxy studied). 

For the wavelengths we chose to use the Johnson- 
Kron- Cousins UBVRI system in the visible and the 
Bessell-Brett system in the near-infrared. Definitions 
of the UBV passbands are given in Buser & Kurucz 
(1978), for R and I in Bcsscll (1990) and for JHK in 
Bessell & Brett (1988). 

To obtain F\(m, Tq) from the theoretical luminosi- 
ties, effective temperatures and gravities along the 
isochrones for each metallicity we preferred to use, 
whenever possible, empirical relations from observed 
stars, rather than theoretical spectra. We performed 
this using bolometric corrections (BC) and T^y-color 
and color-color relations. This was possible for the 
isochrones with metallicities not very different from 
solar: Z=0.008, Z=0.02 and Z=0.05. For lower metal- 
licities, due to the lack of calibration stars, we had 



to take the computed colors of BBCFN, obtained by 
convolving the SEDs of the library of stellar spectra 
of Kurucz with the theoretical passbands. 

Dust analysis. To be able to investigate the pres- 
ence of dust, we have implemented in the code the 
Galactic reddening law of Rieke & Lebofsky (1985) 
which allows us to correct for small amounts of red- 
dening. 

2.3.1. Bolometric corrections 

For early-type stars we calculated the BCs follow- 
ing Code et al. (1976). For late-type stars we used the 
formula of Bessell & Wood (1984) which relates the 
BC to the V — K color. Finally, we used the metal- 
dependent formulae given in Alonso et al. (1995) to 
determine the BCs of the low-MS in the range of tem- 
peratures: 4000-8000K. Because the resulting BCs, 
from different sources, do not match exactly, we have 
interpolated linearly between them to avoid disconti- 
nuities. Fig. 2 summarizes the BC for solar composi- 
tion. 

2.3.2. Temperature- color and color-color conversions 

Where possible we have used empirical studies for 
which the effective temperatures had been determined 
using reliable stellar angular diameter measurements 
(using the lunar occultation technique). For color- 
color conversions we preferred studies which include 
a full set of colors rather than those which give a more 
restricted set, and require the inclusion of subsidiary 
information to complete them. In this way we avoid, 
as far as possible, problems derived from the use of 
different sets of filters. Nowadays the various obser- 
vatories are quite well standardized. It is rather a 
question of different filter sets which were used at dif- 
ferent times. Separate conversions were performed for 
giants and dwarfs: 

• Dwarfs in the range 4000-8000 K: except for U- 
V (for which we used the color-color relations of 
Johnson 1966), we inverted the polynomial fit- 
ting functions of Alonso et al. (1996) to obtain 
each color as a function of the T e ff and the 
metallicity. To obtain their relations they ob- 
served a large sample of MS stars (~ 500) with 
spectral types from F0 to K5 covering a wide 
range in metallicity (-3.0 < [M/H] < 0.5). To 
transform from their J and K filters to our sys- 
tem we used the relations given in Alonso et al. 
(1994). 

• Stars hotter than 8000 K: we first obtained the 
B-Vhom the T eff -(B-V) table of Code et al. 
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(1976) and then used the color-color relations of 
Johnson (1966). We took a smoothly interpo- 
lated transition from Alonso et al. (1996) for 
the boundary at 8000 K. This method works 
since line blanketing for these hot stars is not 
very important. For the case of stars hot- 
ter than 34000 K we used the colors given in 
BBCFN (up to 50000 K obtained by convolv- 
ing with the Kurucz stellar spectra library and 
above 50000 K by assigning pure black-body 
(BB) spectra, see BBCFN for details). Finally, 
to convert from the R and I broad-band filters 
in the Johnson (1966) system to those in the 
Johnson, Kron, Cousins system (Bessell 1990) 
we used the linear relations of Bessell (1983). 

• Dwarfs cooler than 4000 K: due to lack of reli- 
able observations we took the BB effective tem- 
peratures of Johnson (1966). The same paper, 
which does not list any metallicity dependence, 
was used to obtain the different colors by shift- 
ing them to match the well-calibrated metal- 
dependent colors of Alonso et al. (1996) at 
4000 K. In this way we have in fact obtained 
the dependence on metallicity, which becomes 
more important as the temperature decreases, 
as can be seen in Alonso et al. (1996). 

• Giants in the range 3350-4930 K: using the em- 
pirical calibration of Ridgway et al. (1980) we 
first obtained the (V — K) color from T e ff, after 
which we applied the color-color conversions for 
giants of Bessell & Brett (1988) to obtain the 
other colors. For the cooler giants we also used 
the stellar library of Fluks et al. (1994). 

• Giants cooler than 3350 K: unfortunately there 
is a lack of observations for these stars. Espe- 
cially for old stellar populations, the contribu- 
tion from giants with spectral types later than 
M5 is not negligible in red filters especially for 
steep IMFs. For these stars we used the metal- 
dependent model colors of Bessell et al. (1989, 
1991). Since these models do not contain the 
U and B filters, we obtained the U — V and 
B — V colors from the observational stellar li- 
brary of Fluks et al. (1994), taking stars with 
the corresponding V — K. Finally, to avoid prob- 
lems in matching these stars with hotter giants 
we have interpolated the colors between 3200 
and 3350 K. We strongly encourage observers 
to carry out programs to obtain good empirical 
data for cool giants as well as for cool dwarfs. 
It is also important to study empirically the de- 
pendence on the metal content. In this respect 
observations in Baade's window such as those 



obtained by Terndrup et al. (1990 and 1991) 
are very useful. 

The low-mass metal-poor stars. For stars with 
masses lower than 0.6M Q and low metallicities (Z= 
0.0004, Z=0.001 and Z=0.004), we obtained their pa- 
rameters by extrapolating linearly from the low-MS 
range of the isochrones for which we took the colors 
given by the BBCFN. We were forced to this choice 
due to the lack of empirical data to transform from 
the theoretical to the observational plane. 

2.4. Line index synthesis 

As well as the colors, we have also calculated ab- 
sorption line strengths. To make it possible to com- 
pare observations with the line strengths from the 
models, we have worked in the extended Lick-system 
(Worthey et al. 1994). This has the additional ad- 
vantage that many observed spectra are available for 
calibration: stars in that paper and also in Faber et al. 
(1985) and Gorgas et al. (1993). We used the fitting 
functions of Worthey et al. (1994) to relate each line 
index to its corresponding three atmospheric param- 
eters: [M/H], logg and T eff (or V - K). Since for 
the CN1 and CN2 indices these authors do not give 
fitting functions for metallicities below [M/H] = — 1 
we have obtained second order polynomial relations, 
which are given in Table 1. For this we used around 
~ 30 stars of their sample with low metallicities and 
for which they measured the two indices. We also 
have included some features in the near-IR: the Can 
triplet and the Mgi line. We have used the stellar 
spectra of Diaz et al. (1989) to calibrate second or- 
der polynomial relations between each line index and 
its corresponding atmospheric parameters. All these 
fitting functions have been tabulated in Table 1. As 
pointed out by Diaz et al. (1989) it is useful to distin- 
guish between two ranges, metal-poor and metal-rich 
(separated at [M/H] = —0.3). In the first the index 
is strongly dependent on [M/H], while in the second 
it depends much more strongly on gravity. 

We assigned a luminosity class to each live star by 
adopting a rather broad set of criteria to distinguish 
between dwarfs and giants. One of the best discrim- 
inators for this is clearly logg and we have used as 
limits: 



log g > 4.0 
3.5 < log g < 4.0 
log g < 3.5 



dwarfs 
(V-K) < 21og.g-6 
(V-K) >21og ff -6 

giants 



dwarfs 
giants 



(9) 



Because in the gravity range 3.5 < log g < 4, the 
value of log g alone cannot discriminate between the 
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two classes of stars, we decided to use the V — K crite- 
rion cited to divide the subgiant region into two parts. 
Finally, stars with (V- if) < -1.0 or log T eff > 4.63 
are not classified, because they are not important 
for line-strength computations. This equation was 
mainly used when calculating the contributions of the 
M stars to the integrated absorption features. 

The code yields a stellar distribution for a given 
age of the galaxy and metallicity (Z), converted to 
[M/H] by the relation 

[M/H] = log(Z/Z ) (10) 

To calculate an integrated line index the code inte- 
grates the contributions of all the stars (obtained us- 
ing Eq. 11) to calculate the total line index for the 
resulting SSP. An integrated line-index Wi is a flux- 
weighted mean flux and could be expressed either as 
an equivalent width (EW) or in magnitudes. We have 
decided to maintain the convention of Burstein et al. 
(1984), used also in Worthey et al. (1994) of express- 
ing the molecular-band features in magnitudes and 
the atomic-line features in equivalent width (EW) in 
A. The NIR features are expressed in EW as well. We 
calculate Wi in the following way: 

w _ O W{m,T G )N{m,T G )F c {m,T G )dm 

C G N(m,T G )F c (m,T G )dm ( > 

where F c (m,T G ) is the flux in the continuum cor- 
responding to the central wavelength of the spectral 
feature of the star of mass m at T G , obtained by a 
linear interpolation between the closest broad-band 
filter fluxes Fx, and F\ u 

F c (m, T G ) = F Xl + (F\ u F Xl ) (^^) (12) 

where Ai, A; and X u are the central wavelengths of the 
index and the two broad-band filters respectively. 

Finally, to convert the indices expressed in EW to 
magnitudes we use the following relation 

mag = -2.5 log (^=^) (13) 

where AA^ is the wavelength width (in A) of the fea- 
ture. 

2.5. The model input 

To model a single stellar population the main input 
parameters which must be defined are: 

• The IMF shape (unimodal or bimodal as defined 
in Section 2.2), the IMF slope \i and the lower 
and upper cutoff limits m;, m u (we kept them 
fixed at 0.0992 and 72 M ). 



• The age T G of the galaxy (or of the observed 
zone of the galaxy). 

• The metallicity Z . 

2.6. The derived colors and absorption lines 
for a single stellar population 

In Tables 2 and 3 we present the resulting colors 
and line indices for single stellar populations for uni- 
modal and bimodal IMFs with two slopes /x: 1.35 and 
2.35. For each /j, we present the results for different 
metallicities and ages. The differences between the 
two IMFs mainly affect the reddest colors or indices. 
These Tables, as well as others which also include data 
for models with other \x values can be obtained from 
the authors. 
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Index 


a 2 


Z 1 


T tfS 


9 


Z 


T eff 


Zg 


gTeff 


ZT eff 


Constant 


Validity 




CNl 
CNl 


-0.0428 
-0.0089 


0.4230 
0.0818 


-2.6264 
-0.0558 


-7.6603 
-2.8505 


16.2199 
0.0702 


4.5248 
-2.9541 


-0.1189 
-0.0002 


2.1176 
0.7670 


-3.9428 
0.0945 


19.9478 
12.1739 


3980 <T eff 
5100 < T ef/ 


< 5100 

< 11100 


CN2 
CN2 


-0.0196 
-0.0148 


-0.1867 
0.0170 


0.1706 
0.1875 


-0.0520 
-2.6252 


-10.5277 
-0.4784 


-0.1924 
-4.2434 


0.0541 
-0.0485 


0.0734 
0.7027 


2.6526 
0.2019 


-2.5165 
13.5325 


3980 <T ef/ 
5100 < T ef/ 


< 5100 

< 11100 


Caj/l 
Caul 


-0.0065 
0.1269 


0.0024 
-0.8879 


0.4254 
0.6840 


5.5058 
4.9417 


-6.2485 
4.2020 


1.1452 
0.5141 


-0.1162 
-0.3062 


-1.5381 
-1.6106 


1.8636 
-0.7240 


-7.9084 
-7.7763 


Z < -0.3 
Z > -0.3 




Caj/2 
Ca J7 2 


0.1669 
0.2208 


-0.1111 
-2.3469 


0.6564 
0.0913 


9.9490 
1.1574 


-8.8215 
-9.3570 


3.8129 
1.0685 


-0.1269 
-0.3788 


-3.0427 
-0.7946 


2.6914 
3.2014 


-17.0107 
1.9801 


Z < -0.3 
Z > -0.3 




Ca//3 
Ca//3 


0.1335 
0.1801 


-0.1789 
-1.6372 


1.6471 
0.6482 


9.8307 
-1.1326 


-2.5977 
-20.1575 


-4.5822 
-3.4841 


-0.0627 
-0.6283 


-2.9351 
-0.1215 


0.8503 
6.2543 


-0.6660 
10.3075 


Z < -0.3 
Z > -0.3 




Mg/ 
M gI 


0.0425 
0.0302 


0.0453 
0.2442 


0.5193 
-0.8184 


3.9973 
0.9179 


1.9597 
1.3805 


-3.2429 
4.0668 


-0.0278 
0.0241 


-1.1445 
-0.2869 


-0.4173 
-0.3096 


5.9036 
-3.1040 


Z < -0.3 
Z > -0.3 





Table 1: The fitting functions for the CNl and CN2 indices (as defined in Worthey et al. (1994) for [M/H] < -1 
and for the near-IR features as defined in Diaz et al. (1989) (who define the region of the continuum as well). 
Given are the polynomial coefficients where Z should be read as [M/H], g as log g and T e // as logT e ff 
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Z " 0.008 " 0.02 " < 0.0. 



Age 


1 


4 




12 




17 


1 


4 






12 


17 


1 


4 






12 


17 


(M/L) v 


0.75 


2.65 


3.87 


5.37 




7.0 


0.91 


2.96 


5.39 


7.05 


9.14 


1.16 


3.78 


7.20 


10.3 


12.8 


U-V 


0.687 


1.087 


1.302 


1.401 


1 


446 


0.839 


1.328 


1.474 


1 


647 


1.777 


0.955 


1.611 


1.778 


1 


913 


2.063 


B-V 


0.476 


0.780 


0.879 


0.917 





933 


0.603 


0.891 


0.946 


1 


007 


1.049 


0.717 


1.006 


1.058 


1 


100 


1.146 


V-R 


0.354 


0.491 


0.538 


0.557 





574 


0.389 


0.540 


0.569 





600 


0.623 


0.433 


0.600 


0.631 





661 


0.687 


V-I 


0.820 


1.045 


1.140 


1.175 


1 


210 


0.835 


1.144 


1.197 


1 


257 


1.304 


0.867 


1.247 


1.304 


1 


361 


1.411 


V-.I 


1.716 


1.993 


2.136 


2.150 


2 


161 


1.906 


2.287 


2.343 


2 


416 


2.440 


1.997 


2.625 


2.584 


2 


622 


2.640 


V-K 


2.546 


2.843 


3.003 


2.996 


2 


988 


2.815 


3.214 


3.271 


3 


348 


3.356 


2.966 


3.669 


3.580 


3 


597 


3.585 


CN1 


-0.143 


-0.047 


-0.017 


-0.009 


-0 


009 


-0.090 


0.013 


0.026 





043 


0.057 


-0.027 


0.081 


0.101 





117 


0.144 


CN2 


-0.081 


-0.010 


0.016 


0.021 





019 


-0.037 


0.047 


0.056 





074 


0.089 


0.017 


0.113 


0.133 





151 


0.182 


Ca4227 


0.380 


0.950 


1.183 


1.328 


1 


433 


0.601 


1.278 


1.518 


1 


721 


1.884 


0.901 


1.770 


2.021 


2 


234 


2.431 


G-band 


0.698 


3.748 


4.676 


5.037 


5 


317 


1.937 


4.817 


5.398 


5 


742 


5.941 


3.341 


5.673 


6.057 


6 


179 


6.181 


Fe4383 


0.610 


2.958 


3.843 


4.215 


4 


531 


1.951 


4.770 


5.365 


5 


894 


6.300 


3.594 


6.605 


7.236 


7 


614 


8.004 


Ca4455 


0.519 


1.066 


1.276 


1.368 


1 


442 


0.907 


1.515 


1.655 


1 


795 


1.898 


1.382 


2.005 


2.134 


2 


233 


2.332 


Fe4531 


1.594 


2.596 


2.940 


3.096 


3 


228 


2.214 


3.244 


3.482 


3 


715 


3.898 


2.877 


3.915 


4.158 


4 


339 


4.531 


Fc4668 


0.674 


2.237 


2.817 


2.950 


2 


969 


2.532 


4.640 


4.910 


5 


293 


5.547 


5.029 


7.608 


8.008 


8 


315 


8.702 


H/3 


4.909 


2.390 


1.853 


1.613 


1 


402 


4.156 


2.007 


1.625 


1 


357 


1.157 


3.267 


1.656 


1.339 


1 


137 


0.982 


Fg5015 


2.948 


4.105 


4.521 


4.645 


4 


652 


3.965 


5.301 


5.440 


5 


645 


5.767 


5.117 


6.443 


6.539 


6 


654 


6.772 


M gl 


0.019 


0.054 


0.071 


0.081 





092 


0.034 


0.089 


0.106 





123 


0.136 


0.064 


0.136 


0.155 





171 


0.189 


Mg2 


0.084 


0.151 


0.184 


0.202 





216 


0.123 


0.218 


0.247 





274 


0.294 


0.172 


0.295 


0.327 





353 


0.377 


Mg6 


1.560 


2.474 


2.958 


3.192 


3 


333 


2.056 


3.391 


3.790 


4 


052 


4.233 


2.642 


4.316 


4.703 


4 


972 


5.123 


Fe5270 


1.338 


2.166 


2.448 


2.573 


2 


673 


1.975 


2.827 


2.995 


3 


176 


3.321 


2.613 


3.447 


3.623 


3 


760 


3.918 


Fe5335 


1.004 


1.794 


2.062 


2.181 


2 


293 


1.670 


2.530 


2.705 


2 


884 


3.016 


2.483 


3.377 


3.547 


3 


689 


3.830 


Fe5406 


0.573 


1.121 


1.317 


1.395 


1 


471 


1.001 


1.624 


1.746 


1 


885 


1.987 


1.508 


2.183 


2.304 


2 


408 


2.524 


Fe5709 


0.537 


0.726 


0.785 


0.794 





798 


0.751 


0.970 


0.984 


1 


014 


1.022 


1.018 


1.219 


1.229 


1 


226 


1.243 


Fe5782 


0.356 


0.589 


0.653 


0.666 





672 


0.583 


0.818 


0.836 





882 


0.908 


0.807 


1.046 


1.064 


1 


083 


1.130 


NaD 


1.306 


2.119 


2.438 


2.703 


2 


979 


1.941 


2.994 


3.402 


3 


747 


4.093 


2.835 


4.194 


4.678 


5 


093 


5.496 


TiOj 


0.028 


0.034 


0.038 


0.039 





041 


0.026 


0.040 


0.042 





045 


0.048 


0.024 


0.049 


0.053 





057 


0.060 


TiO„ 


0.028 


0.047 


0.055 


0.058 





062 


0.025 


0.061 


0.067 





074 


0.080 


0.027 


0.089 


0.096 





105 


0.112 


Caj;l 


1.212 


1.324 


1.364 


1.383 


1 


406 


1.731 


1.781 


1.772 


1 


786 


1.793 


1.813 


1.961 


1.901 


1 


887 


1.892 


Ca J; 2 


3.389 


3.589 


3.672 


3.714 


3 


759 


4.454 


4.576 


4.502 


4 


502 


4.455 


4.542 


4.815 


4.592 


4 


496 


4.439 


Ca„3 


2.850 


3.037 


3.114 


3.162 


3 


213 


3.865 


3.868 


3.745 


3 


714 


3.631 


3.955 


4.088 


3.752 


3 


582 


3.477 


Mg„ 


0.491 


0.634 










0.641 












0.843 


1.018 




1 




1.076 











Z 








0.008 














0.02 
















0.05 










Age 


1 




4 


8 




12 


17 




1 


4 


8 




12 


17 




1 




4 


8 




12 




17 


(M/L) v 


4.82 




12.0 


15.5 




19.6 


23.6 




3.06 


14.3 


22.3 


27.0 


32.5 


7.59 




18.3 


30.0 


39.7 


46.2 


U-V 


0.668 


1 


114 


1.346 


1 


465 


1.530 





807 


1.334 


1.518 


1 


701 


1.842 





920 


1 


599 


1.801 


1 


951 


2 


103 


B-V 


0.470 





800 


0.905 





953 


0.978 





593 


0.900 


0.975 


1 


041 


1.089 





703 


1 


007 


1.076 


1 


127 


1 


177 


V-R 


0.358 





521 


0.574 





600 


0.624 





397 


0.568 


0.615 





651 


0.679 





440 





624 


0.674 





714 





744 


V-I 


0.832 


1 


135 


1.244 


1 


302 


1.353 





872 


1.235 


1.340 


1 


416 


1.480 





915 


1 


337 


1.453 


1 


544 


1 


607 


V-J 


1.693 


2 


113 


2.275 


2 


330 


2.377 


1 


886 


2.382 


2.513 


2 


608 


2.668 


2 


007 


2 


702 


2.772 


2 


867 


2 


921 


V-K 


2.492 


2 


958 


3.139 


3 


184 


3.222 


2 


756 


3.295 


3.430 


3 


532 


3.582 


2 


927 


3 


714 


3.744 


3 


825 


3 


860 


CN1 


-0.151 


-0 


060 


-0.034 


-0 


031 


-0.033 


-0 


101 


-0.002 


0.010 





023 


0.033 


-0 


040 





068 


0.088 





101 





121 


CN2 


-0.089 


-0 


024 


-0.003 


-0 


002 


-0.006 


-0 


048 


0.031 


0.039 





052 


0.063 





002 





098 


0.119 





133 





160 


Ca4227 


0.457 


1 


139 


1.424 


1 


624 


1.774 





665 


1.415 


1.735 


1 


977 


2.182 





927 


1 


826 


2.148 


2 


405 


2 


620 


G-band 


0.549 


3 


681 


4.557 


4 


880 


5.115 


1 


840 


4.716 


5.256 


5 


546 


5.691 


3 


303 


5 


603 


5.928 


6 


004 


5 


974 


Fe4383 


0.509 


2 


971 


3.836 


4 


201 


4.494 


1 


817 


4.684 


5.293 


5 


780 


6.136 


3 


461 


6 


488 


7.118 


7 


450 


7 


786 


Ca4455 


0.511 


1 


117 


1.343 


1 


455 


1.543 





890 


1.527 


1.698 


1 


843 


1.952 


1 


351 


1 


986 


2.132 


2 


235 


2 


330 


Fc4531 


1.603 


2 


709 


3.084 


3 


278 


3.435 


2 


200 


3.290 


3.583 


3 


826 


4.018 


2 


830 


3 


899 


4.177 


4 


362 


4 


545 


Fe4668 


0.454 


1 


993 


2.532 


2 


639 


2.648 


2 


250 


4.365 


4.651 


5 


001 


5.222 


1 


689 


7 


316 


7.753 


8 


030 


8 


383 


H/3 


4.891 


2 


169 


1.601 


1 


318 


1.081 


4 


133 


1.860 


1.401 


1 


114 


0.894 


3 


267 


1 


579 


1.209 





983 





824 


Fe5015 


2.771 


3 


977 


4.375 


4 


488 


4.497 


3 


793 


5.129 


5.277 


5 


462 


5.565 


4 


956 


6 


274 


6.378 


6 


476 


6 


580 


Mgl 


0.024 





068 


0.088 





101 


0.115 





037 


0.098 


0.121 





139 


0.154 





061 





138 


0.162 





180 





198 


Mg 2 


0.090 





167 


0.203 





224 


0.240 





125 


0.228 


0.264 





292 


0.313 





168 





298 


0.337 





364 





388 


Mgb 


1.575 


2 


581 


3.042 


3 


262 


3.387 


2 


056 


3.444 


3.857 


4 


100 


4.262 


2 


631 


4 


342 


4.749 


5 


008 




147 


Fe5270 


1.324 


2 


236 


2.527 


2 


667 


2.773 


1 


928 


2.827 


3.024 


3 


198 


3.332 


2 


529 


3 


401 


3.595 


3 


722 


3 


861 


Fe5335 


1.015 


1 


883 


2.162 


2 


298 


2.415 


1 


635 


2.529 


2.725 


2 


894 


3.013 


2 


401 


3 


313 


3.491 


3 


620 


3 


739 


Fe5406 


0.572 


1 


175 


1.379 


1 


472 


1.553 





962 


1.612 


1.752 


1 


884 


1.976 


1 


434 


2 


127 


2.255 


2 


349 


2 


448 


Fe5709 


0.480 





642 


0.683 





673 


0.659 





684 


0.869 


0.849 





858 


0.843 





943 


1 


117 


1.089 


1 


055 


1 


051 


Fe5782 


0.337 





562 


0.623 





633 


0.637 





541 


0.769 


0.784 





820 


0.837 





743 





981 


0.993 


1 


000 


1 


035 


NaD 


1.690 


2 


895 


3.334 


3 


743 


4.121 


2 


253 


3.633 


4.301 


4 


727 


5.153 


2 


981 


4 


608 


5.320 


5 


842 


6 


262 


TiOj 


0.029 





041 


0.047 





051 


0.055 





030 


0.049 


0.056 





062 


0.068 





029 





057 


0.067 





075 





080 


TiOjj 


0.031 





061 


0.074 





082 


0.090 





034 


0.079 


0.096 





108 


0.120 





036 





104 


0.124 





140 





152 


C« ; ,l 


1.220 


1 


377 


1.424 


1 


454 


1.483 


1 


718 


1.807 


1.823 


1 


844 


1.859 


1 


767 


1 


936 


1.906 


1 


906 


1 


917 


Ca„2 


3.451 


3 


780 


3.881 


3 


957 


4.024 


4 


286 


4.451 


4.392 


4 


397 


4.366 


4 


324 


4 


593 


4.383 


4 


296 


4 


254 


Ca 7J 3 


2.936 


3 


246 


3.341 


3 


422 


3.494 


3 


606 


3.618 


3.479 


3 


448 


3.375 


3 


580 


3 


674 


3.298 


3 


117 


3 


022 


Mgj-j 


0.574 





786 


0.850 





909 


0.960 





691 


0.876 


0.944 





978 


1.007 





889 


1 


082 


1.134 


1 


170 


1 


189 



Table 2: The model observables for SSPs with a unimodal IMF with slopes 1.35 and 2.35. The age is in Gyr. 
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1.35 


1.35 


1.35 


Z 


















0.02 


















Age 


1 




4 




12 


17 


1 


4 




12 


17 


1 


4 




12 




17 


(M/L) v 


0.48 




1.60 


2.27 


3.10 


3.97 


0.59 


1.80 


3.19 


4.11 
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Table 3: The model observables for SSPs with a bimodal IMF with slopes 1.35 and 2.35 
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Fig. 3. — Comparison of the colors predicted in our SSP models with those obtained by different authors versus 
age and metallicity for a Salpctcr IMF. 
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2.6.1. Comparison with other authors 



As a check on our method we compared our syn- 
thetic colors with those published by Buzzoni (1989) 
(we selected those models which include the red hor- 
izontal branch and have a mass-loss rate parameter 
r] = 0.3), Bruzual & Chariot (1996) (private commu- 
nication, a revised version of the results of Bruzual 
& Chariot 1993), Tantalo et al. (1996) (which in- 
cludes the revised results of Bressan et al. 1994) and 
W94. In Fig. 3 we present the color-age (for solar 
metallicity) and color-metallicity diagrams for young 
(2 Gyr) and old (12 Gyr) stellar populations. In most 
of the plots we see that our results are always in better 
agreement with the values obtained by Tantalo et al. 
(1996) and Bruzual & Chariot (1996), while we differ 
substantially from the results obtained by W94 and 
Buzzoni (1989). We also note that the older mod- 
els of the Padova group (Bressan et al. 1994) give 
much redder V — K colors than we find (similar to 
colors of W94) while the model of Bruzual & Char- 
lot (1993) is much bluer than ours for moderate ages. 
The similarity between our results and those of Tan- 
talo et al. (1996) is due to the use of the same the- 
oretical isochrones, whereas the small differences are 
caused by the different conversions to the observa- 
tional plane. The very red V — K colors of W94 for 
cool stars are attributed to his bolometric corrections. 
A nice paper analyzing the differences between these 
models has been written by Chariot et al. (1996). 
They found that the scatter between authors is due 
to the use of different isochrones rather than spectral 
calibrations. Finally, the colors of Buzzoni (1989) dif- 
fer very much from the others. 

For the absorption lines we have compared our re- 
sults with those of W94, Bruzual & Chariot (1996) 
and Bressan et al. (1996). In Figs. 4 and 5 we in- 
vestigate the behavior of representative indices, from 
eight different elements. In Fig. 4 these indices are 
plotted as a function of age for solar composition. 
In Fig. 5 we plot the same synthetic indices versus 
the metallicity for young (2 Gyr) and old (12 Gyr) 
stellar populations. Looking at these figures we infer 
that the agreement is generally good. Once again our 
models agree better with Bruzual & Chariot (1996) 
than with W94. In particular we can notice that our 
line-strength predictions are in general stronger than 
the ones of W94, except for H/3, as expected, and for 
TiO/. Since most of these indices are calculated di- 
rectly from the theoretical isochrone parameters with- 
out any intermediate conversion to the observational 
plane, as explained in Section 2.4, we attribute this 
to the fact that our isochrones are slightly cooler, 
especially for the MS and the turnoff stars, as well 



5 10 15 5 10 15 




5 10 15 5 10 15 

Age (Gyr) Age (Gyr) 



Fig. 4. — A representative set of our synthetic indices 
versus age for a Salpeter IMF and solar composition 
compared with the same indices calculated by differ- 
ent authors (from Bressan et al. (1996) only H/3, Mg2 
and < Fe >= Fe5270 + Fe5335 are available). 
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Fig. 5. — The synthetic indices versus metallicity for a 
Salpeter IMF and for assumed ages of 2 and 12 Gyr, 
compared with the equivalent plots from W94 and 
Bressan et al. (1996). 
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as the upper- RGB and AGB stars. This difference 
in the temperature is mainly caused by the adopted 
opacities. However other different stellar evolution 
prescriptions are also significant (see for more details 
Chariot et al. 1996). The disagreement among the 
models is greatest for NaD, Ca4227 and TiO/. This 
is well understood if we look at the index- (V — K) 
diagrams of Gorgas et al. (1993) and Worthey et al. 
(1994) index-(y — K) plots for observed stars. The 
dependence of Ca4227 and NaD on temperature is 
very steep (compared with e.g iron indices) for EWs 
of ~ 1 and ~ 2 respectively (where differences in the 
integrated predictions also become appreciable), es- 
pecially for dwarf stars. The case of TiO/ is slightly 
more difficult to understand, since for very cool stars 
this index changes rapidly as a function of tempera- 
ture for both, dwarfs and giants. Worthey 's TiO/ is 
stronger than ours because his bolometric corrections, 
the same reason that his V — K is redder. This is clear 
when one realises that both are strongly dependent on 
very red stars, for which W94 has a problem with the 
BCs (see Chariot et al. 1996). 

2.6.2. Comparison with observed clusters 

Since clusters are considered to be coeval, they are 
ideal for testing the synthetic colors and absorption 
lines of our models. In Fig. 6 we have compared our 
model colors and absolute magnitudes with the Mr 
vs. V — K color-magnitude diagrams of Houdashclt 
et al. (1992) for M 67 and Frogel et al. (1981) for 
47 Tuc. For M 67 we used an isochrone of 4.5 Gyr and 
solar metallicity, while for 47 Tuc we used an age of 
17 Gyr and Z=0.004. Notice that for 47 Tuc our colors 
are the same as those of BBCFN, since there is no 
difference in the calibration for low metallicities. The 
fit for M 67 is also excellent, in general the agreement 
in V — K is better than 0.05 mag everywhere. W94 
obtained a good fit for Z=0.0134 and 4.5 Gyr, while 
we use Z=0.02. This difference reflects the difference 
in temperature of the isochrones in general. 

Next, we made a comparison between our models 
and the integrated colors as well as line-strengths of 
some globular clusters. Two color-color diagrams and 
V — K vs. Mg2 are given in Fig. 7, while in Fig. 8 a 
number of index-Mg 2 diagrams are represented. For 
the colors, the observational data of Galactic and 
M31 globular clusters are from Burstein et al. (1984), 
who used raw data from different authors (references 
therein). We selected, following the classification of 
Searle et al. (1980), UBV data for the LMC glob- 
ular clusters from Bica et al. (1992) and Girardi et 
al. (1995). The line-strengths of Galactic and M 31 
globular clusters are from Burstein et al. (1984) and 
for other Galactic globular clusters from Covino et 



al. (1995). The integrated model observables have 
been obtained with a bimodal IMF of slope \i = 1.35 
for SSPs with different metallicities and ages (ranging 
from 1 to 17 Gyr). In Fig. 8 we see that our models fit 
the data very well and that the required metallicity is 
always lower than solar. As one can see in Fig. 8 our 
models also fit very well the plots of G band, H/3 and 
iron indices vs. Mg 2 . For instance, the H/3-Mg 2 di- 
agram indicates that the ages of these clusters must 
be older than say ~ 10 Gyr. On the other hand, 
the NaD and CN1 features do not fit so well. The 
NaD mismatch can be probably attributed to inter- 
stellar absorption, which can increase the EW of the 
observed features by up to 1 A (Gorgas et al. 1993). 
However, for the most metal-rich clusters, we some- 
times observe differences even higher than 1 A. Bica 
et al. (1991) found similar results for a large sample 
of galaxies. For the CN lines the mismatch cannot be 
attributed to extinction since the line is purely stellar. 
Here only for the most metal-rich galaxies is the mis- 
match large. There are all clusters of M31. It could 
be that abundance ratio variation effects play a role 
here in the metal-rich galaxy M 31 (W94). W94 also 
attributed the CN mismatch observed in the metal- 
rich galaxy M 31 to abundance ratios effects. This 
effect could also be responsible for part of the NaD 
mismatch. 
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Fig. 6. — Color-magnitude diagram of M 67 and 
47 Tuc together with our stellar colors. We plot stars 
of M 67 from Houdashelt et al. (1992) and stars of 
47 Tuc by Frogcl et al. (1981). 



Fig. 8. — Comparison of our integrated indices with 
observations of Galactic and M31 globular clusters in 
a number of index-Mg2 diagrams. Our indices are 
obtained with a bimodal IMF of slope /i = 1.35 for 
SSPs of different metallicities, ranging in age from 1 
to 17 Gyr. 



Fig. 7. — Comparison of our models with optical and 
IR colors and Mg 2 for Galactic, M31 and LMC globu- 
lar clusters (see the text for more details). Presented 
are models with a bimodal IMF of slope \i= 1.35 for 
SSPs of different metallicities, with age ranging from 
1 to 17 Gyr. 
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3.4 



3.2 



WithouL HB /AGE 
Including all 
Shift +50K 



9 



16 Cyr 



8 Gyr 



Fig. 9. — V — K for models with and without 
HB+AGB, and shifting the isochrone by 50 K. 



! .2 



O Without HB/AGB 
■ Including all 
x Shift +50K 



16 Gyr 



8 Gyr 



0.02 



0.04 



0.06 



Fig. 10. — U — V for models with and without 
HB+AGB, and shifting the isochrone by 50 K. 
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2.7. Errors in integrated line indices and col- 
ors 

Since models of this kind are extremely compli- 
cated, the determination of the errors in colors and 
lines is very difficult. Part of the physics involved is 
not completely understood (e.g. molecular opacities), 
and even where it is understood calculations can be 
very difficult (e.g. accounting for convection). Er- 
rors can be due to uncertainties in theoretical stellar 
parameters, especially for advanced stages of stellar 
evolution, to the slope of the IMF, about which little 
theoretical is known, to the conversion from the theo- 
retical to the observational plane using observed and 
synthetic stars, and to many more factors. 

One way to estimate our errors is to compare our 
results with results from other authors (see Section 
2.6.1). But since some errors are systematic, and 
possibly different for each group, this is not enough. 
For this reason we also present the differences in in- 
tegrated colors and indices between the values for 
whole model, and those of all stages except the HB 
and the AGB for which the isochrone calculations 
might contain large uncertainties. In Fig. 9 we present 

V — K with and without HB+AGB for various ages 
and metallicities. Without the latest phases V — K is 
quite well-behaved, i.e. monotonic in age and metal- 
licity. In Fig. 10, we see that the advanced stages 
contribute considerably in U — V, which could maybe 
explain why our models are amongst the reddest in 
this color. 

The differences are quantified in Table 4. In the 
last column we give the maximum differences, in mag- 
nitude, or as a fraction of the line strength. These 
numbers have the same order of magnitude as the 
differences between the authors in Fig. 3. 

This table gives an idea of the relative accuracy of 
the various colors and lines. The possibilities for sys- 
tematic errors are considerable. The worst colors are 
U-V and V-K, while the models for lines CNi, CN 2 , 
Ca4227, Fe4668 and Fe5782 are the most uncertain. 
These possible errors have to be taken into account 
when applying the model to observational data. Note 
that these are not (and do not have to be) the lines 
for which the comparison between various authors is 
the worst. 

As an additional test we have shifted the whole 
isochrone to the blue by 50 K, which might be a rea- 
sonable uncertainty in the isochrones. The results are 
also plotted in Figs. 9 and 10. As expected, the ef- 
fect is systematic but does not change either U-V or 

V — K by more than 0.05 mag. 

The comparison with M 67 or 47 Tuc shows that 
here our error in V — K is probably smaller than 



0.05 mag everywhere. Assuming that the whole RGB 
(and further stages) are 0.1 mag too blue (or too red) 
we find that the integrated V — K is 0.07 mag too blue 
(or too red) for several ages and metallicities. This 
means that the error due to RGB mismatch is not 
larger than 0.04 mag, although at Z=0.05 the errors 
might be larger, since no globular clusters are present 
here for calibration. 

In order to quantify how any mismatch in the cal- 
ibration of the K vs. V — K diagrams of red gi- 
ant branches of globulars (see Fig. 6) is propagated 
through the models, we assumed that the stars on the 
whole RGB (and further stages) are 0.1 mag too blue 
or 0.1 mag too red. We then found that the integrated 
V — K colors are 0.07 mag too blue with respect to 
the red propagation for models of 17 Gyr, and Z=0.02 
and 0.05. For models of 1 Gyr and the same metallici- 
ties we obtained a difference of 0.08 mag. This means 
that we may assume that the error in integrated V—K 
up to the tip of the RGB is smaller than 0.04 mag, 
although at Z=0.05 this error could be larger, since 
no globular clusters here are present to calibrate the 
models. At low metallicities there are no differences 
between our integrated colors and those of BBCFN. 

2.7.1. The effect of the very low mass stars on the 
integrated observables 

In this section we briefly quantify the effects of very 
low-mass stars (< 0.6Mq) on the integrated observ- 
ables. Since their effects are more important in the 
red we have calculated the effects on V — K and NaD. 
In Table 5 we present integrated V — K and NaD for 
models with a unimodal IMF, a bimodal IMF (which 
diminishes the contributions of these very low-mass 
stars) and an IMF with a lower-mass cutoff of pre- 
cisely 0.6M Q . As one expects their influence is larger 
for steeper IMF slopes and greater ages. See also 
Fig. 14. 
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Color /Index 


mag 


Z=0.02, 8 Gyr 


Z=0.05, 8 Gyr 


Z=0.02, 16 Gyr 


Z=0.05, 16 Gyr 


max. rcl. error 


U - V 


y 


0.156 


0.207 


0.121 


0.162 


0.207 


B - V 


y 


0.060 


0.085 


0.039 


0.060 


0.085 


V - R 


y 


0.019 


0.035 


0.003 


0.014 


0.035 


V - I 


y 


0.030 


0.063 


-0.010 


0.017 


0.063 


V - ,7 


y 


0.115 


0.122 


-0.054 


0.009 


0.122 


V - K 


y 


0.163 


0.170 


-0.070 


0.017 


0.170 


CNi 


y 


0.029 


0.041 


0.033 


0.040 


0.041 


CN 2 


y 


0.031 


0.045 


0.037 


0.043 


0.045 


Ca4227 


n 


0.028 


0.125 


-0.033 


0.058 


0.125 


G4300 


n 


0.043 


0.010 


0.030 


0.004 


0.043 


Fo4383 


n 


0.079 


0.063 


0.045 


0.035 


0.079 


Ca4455 


n 


0.071 


0.079 


0.038 


0.053 


0.079 


Fc4531 


n 


0.050 


0.064 


0.023 


0.036 


0.064 


Fc4668 


n 


0.132 


0.113 


0.103 


0.085 


0.132 


H/3 


n 


-0.067 


-0.088 


0.000 


-0.007 


0.088 


Fc5015 


n 


0.067 


0.068 


0.051 


0.059 


0.068 


Mgi 


y 


0.010 


0.025 


0.000 


0.014 


0.025 


Mg 2 


y 


0.011 


0.027 


-0.003 


0.012 


0.027 


Mg b 


n 


0.009 


0.023 


-0.030 


-0.010 


0.030 


Fo5270 


n 


0.059 


0.070 


0.029 


0.039 


0.070 


Fo5335 


n 


0.056 


0.073 


0.019 


0.041 


0.073 


Fo5406 


n 


0.086 


0.103 


0.042 


0.062 


0.103 


Fe5709 


n 


0.099 


0.100 


0.088 


0.092 


0.100 


Fe5782 


n 


0.125 


0.165 


0.094 


0.126 


0.165 


NaD 


n 


-0.004 


0.045 


-0.051 


-0.013 


0.051 


TiOj 


y 


0.000 


0.005 


-0.002 


0.002 


0.005 


TiOjj 


y 


0.001 


0.015 


-0.004 


0.009 


0.015 


Cajjl 


n 


0.030 


0.065 


0.015 


0.053 


0.065 


Ca /7 2 


n 


0.043 


0.065 


0.014 


0.053 


0.065 


Ca J7 3 


n 


0.062 


0.107 


0.026 


0.098 


0.107 


M gI 


n 


0.009 


0.002 


-0.033 


-0.026 


0.033 



Table 4: Predicted differences in integrated colors and indices for SSPs with and without the inclusion of the 
HB+AGB phases of stellar evolution. The given numbers indicate the relative errors in magnitude, or as a fraction 
of the line-strength 















M = 


0.35 


















V - 


- K 










NaD 










1 Gyr 






16 Gyr 






1 Gyr 






16 Gyr 




z 


U 


B 


O 


U 


B 


O 


U 


B 


O 


U 


B 


O 


0.008 


2.70 


2.70 


2.70 


2.94 


2.94 


2.90 


1.29 


1.28 


1.27 


2.46 


2.44 


2.13 


0.02 


2.97 


2.97 


2.97 


3.30 


3.30 


3.27 


1.97 


1.97 


1.95 


3.56 


3.54 


3.26 


0.05 


3.11 


3.11 


3.11 


3.54 


3.54 


3.50 


2.95 


2.95 


2.93 


5.08 


5.07 


4.83 














M = 


1.35 


















V - 


- K 










NaD 










1 Gyr 






16 Gyr 






1 Gyr 






16 Gyr 




Z 


U 


B 


o 


U 


B 


O 


U 


B 


O 


U 


B 


O 


0.008 


2.55 


2.54 


2.53 


2.99 


2.94 


2.87 


1.31 


1.28 


1.22 


2.92 


2.71 


2.17 


0.02 


2.82 


2.81 


2.80 


3.34 


3.29 


3.23 


1.94 


1.92 


1.86 


4.00 


3.83 


3.33 


0.05 


2.97 


2.96 


2.95 


3.59 


3.54 


3.46 


2.83 


2.82 


2.77 


5.43 


5.32 


4.92 














M = 


2.35 


















V - 


- K 










NaD 










1 Gyr 






16 Gyr 






1 Gyr 






16 Gyr 




Z 


U 


B 


O 


U 


B 


O 


U 


B 


O 


U 


B 


O 


0.008 


2.49 


2.38 


2.33 


3.22 


2.97 


2.84 


1.69 


1.43 


1.21 


4.05 


3.19 


2.23 


0.02 


2.76 


2.65 


2.60 


3.56 


3.32 


3.20 


2.25 


2.04 


1.83 


5.05 


4.31 


3.43 


0.05 


2.93 


2.82 


2.76 


3.86 


3.57 


3.43 


2.98 


2.85 


2.68 


6.19 


5.72 


5.04 



Table 5: Test of the influence of the stars of masses lower than 0.6 M Q on the integrated V — K color and NaD 
index for two ages (1 and 16 Gyr) and three metallicities (Z=0.008, Z=0.02 and Z=0.05). 'U' means that the 
observables were obtained using a unimodal IMF, which includes all the stars up to 0.0992 Mq, 'B' means that we 
used a bimodal IMF which diminishes the influence of these very low-mass stars, and 'O' means that we avoided 
the low-mass stars by taking an IMF with a lower mass-cutoff of 0.6 M Q . As expected the influence of such stars 
is larger for steeper IMF slopes and greater ages 
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3. The chemo-evolutionary model 

In this section we present a new and fully elabo- 
rated spectrophotomctric stellar population synthesis 
model which, apart from an assumed IMF, needs only 
an analytical functional form for the SFR (see 3.1) 
to follow the evolution of a galaxy from an initial gas 
cloud to the present time, including the chemical evo- 
lution by taking into account mass loss from stars and 
supcrnovae. We explain the ingredients of the models 
together with the method of calculation in 3.2, while 
all the details about the ejecta are explained in 3.3. 
The output from the code will be effectively a mix- 
ture of SSPs of different ages and metallicities, which 
have been calculated in the way described in Section 
2. We provide here a discussion of the influence of 
each of the main input parameters and we compare 
our chemical evolution model with the model of Ari- 
moto & Yoshii (1986) (hereafter AY86). Finally we 
present some tables with model observables obtained 
for a set of representative input parameters. 

3.1. The Star Formation Rate 

To quantify how much gas is converted to stars at 
each time, we need to define the Star Formation Rate 
(SFR) C(t). We assume that it is proportional to a 
power, k, of the fractional gas mass f g (t) (f g (t) = 
M g (t)/Mc(t) where M g (t) is the mass of the gas at 
time t and Ma(t) is the total mass of the zone at that 
time, so that M G (t) = M g {t) + M s (t) where M s (t) is 
the stellar mass, including remnants. Then 

C(t) = V f g {t) k (14) 

The power k is observationally estimated at between 
1 and 2 according to Schmidt (1959), or k = 1.3 ±0.3 
(Kennicutt 1989). In this paper we take k = 1. v is a 
constant, which fixes the timescale of star formation 
(see Arimoto & Yoshii, 1986,1987) and should depend 
on the overall physical state of the gas: temperature, 
density, and magnetic field strength etc. 

3.2. The stellar population synthesis and the 
chemical evolution 

We broadly follow the mathematical formalism es- 
tablished by AY86 and Casuso (1991). We assume a 
fixed volume, a well-defined zone, in the galaxy under 
study, initially containing gas, with mass M g (t = 0) 
which in the present paper is taken as fixed, i.e. no 
gas flows into or out of the volume considered. When 
the physical conditions appropriate for gravitational 
collapse are reached stars are formed, with the appro- 
priate SFR and the IMF. 



The implementation is as follows: during a speci- 
fied time interval At, a generation of stars, with prop- 
erties determined by the SFR, the IMF and the metal- 
licity at that moment, is created if f g (t) is greater 
than a certain lower limit f gmin - Successive genera- 
tions of stars will be formed until the present time 
(or a specified epoch Tq) so that the light which is 
observed comes from a composite set of stars from all 
the generations, (SSPs), each with its corresponding 
age and metallicity, surviving at time Tq. 

3.2.1. The stellar population model 

We start by assuming that all the mass is in 
gaseous form, i.e. zero initial mass fraction in stel- 
lar form (M 3 (t = 0) = 0), therefore f g (t = 0) = 1, 
and that the gas has a certain metallicity (usually 
Z(t = 0) = 0). The evolution of the fractional gas 
mass f g (t) is governed by the differential equation 

#M = -C{t)+F(t) (15) 

where F(t) is the fractional mass of gas (expelled by 
stars and due to inflow from outside the zone) enter- 
ing the zone per unit time at time t, during the last 
time-interval (t, t — At). We first must calculate F(t), 
which comprises the total ejecta (from the stars which 
had been created before the given epoch) in the most 
recent time- interval (t, t — At), 

rt ,ro«j(t'-At) 

F(t) = / B(m,t-t')R z (m,t')dmdt'+P(t) 

J At Jm w (t') 

(16) 

where the integration over time represents the contri- 
butions to the fractional gas mass of the stars of each 
SSP created until t. A given SSP has an age of t' at 
time t. Therefore, the last contributing SSP has an 
age t' = At and the oldest one is the first SSP, with 
an age t' = t. In this equation: 

B(m, t — t') (called the birth function) is the gas mass 
fraction per unit time per unit mass, which goes into 
stars of mass m at time t — t' , divided by the total 
mass of the zone. 

R z (m, t') is the gas mass fraction ejected by a star of 
mass m and metallicity Z until the time t', or during 
its lifetime t m if t' >t m . 

m w (t') is the lowest mass corresponding to a post- 
RGB star which is ejecting matter via winds at the 
age t'. 

md(t' — At) is the mass of the star whose lifetime t m 
is equal to the age t' — At. 

P(t) is a term which represents the net inflow of ma- 
terial entering the zone. 
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Splitting the integral over the mass into two terms 
we get: 



F(t) 



+P(t) 



rt rm d (t') 

I I B(m,t-t')R z (m,t')dm+ 

J At Jm„,(t') 



I At Jm m (t') 
r m d (t'-At) 



Jm d (t') 



B(m, t — t')R z (m, t m )dm 



dt' 
(17) 



where the first term represents the contribution to 
the total ejecta from stars in post-RGB stages, whose 
matter is ejected into the interstellar medium via 
winds, and the second term represents the contri- 
bution from the stars which at time t' had already 
reached the end of their evolution and ejected the 
whole of their initial masses except their remnants. 
The masses m w (t'), md(t') and rrid(t' — At) are ob- 
tained from the corresponding isochrones. We note 
that these quantities are approximate, because they 
have been obtained from the closest isochrone to the 
required age which has the closest metallicity to that 
of the desired SSP. However the fractions R z (m,t') 
are obtained directly from the isochrones (see Section 
2.1.2) and the R z (m 7 t m ) from the final ejecta (see 
Section 3.3). 

For analytical simplicity we separate B(m,t) into 
time-dependent and mass-dependent terms, C(t) and 
$(m), yielding 



performed numerically and corrects errors generated 
by the need to discretize the time evolution. 

The finite time step also leads to the following 
problem. Consider a group of post-RGB stars which 
at time t are ejecting matter into the ISM so that their 
contributions are included in F(t). Let us now assume 
that at some point during the next step t + At these 
stars die, so that their total contributions R z (m,t m ) 
are included in F(t + At). Therefore F(t + At) also 
contains the part of the ejected matter previously 
computed. However at this time t + At this gas, could 
for example go into the formation of new stars but at 
the same time was computed as a newly ejected mat- 
ter. To correct it, we subtract this gas from F(t+At). 

Infall and Outflow. In this paper we consider a 
closed-box model in which the zone does not suffer 
any kind of interchange of gas with its neighborhood. 
However, to be more general the equations written 
above include this possibility by means of the term 
P(t) which can have different forms (see for example 
Lacey & Fall (1985), or Clayton (1985,1986). Obvi- 
ously, for a closed-box model the term P(t) is 0. 

3.2.2. Chemical evolution 

Historically, see for example Tinsley (1980) and 
AY86, the total fraction of the heavy elements ejected 
by a star of initial mass m was written as 



B(m,t) = ${m)C(t) (18) 

We have solved the equation (15) without using 
the approximation of instantaneous recycling (Tins- 
ley 1980), but evaluating f g (t) at time t using the pre- 
vious values of f g (t - At), C(t - At) and F(t - At) 
via 

f g (t) = f g (t-At) + ^At (19) 

If the resulting value of f g (t) is greater than a 
pre-set value f gmin the SFR, C(t), is evaluated from 
Eq. 15, but if it is smaller no stars are being formed 
and therefore C(t) is set to 0. We apply this thresh- 
old based on a suggestion by Kennicutt (1989) that 
below a certain critical gas density there is no massive 
star formation at all. However, Caldwell et al. (1994) 
conclude the opposite. For that reason we decided for 
the present to set this free parameter f 9min to 0. 

Mass Conservation. In order to guarantee mass 
conservation in our calculations, we have normalized 
to unity the constant /3 at each time-step At: this is 



E z (m) = Q z (m)+Z g (t-t m )(R z (m,t)-Q z (m)) (20) 

where Z g represents an average metallicity Z or that 
of any individual element Fe, O, etc.; Z g (t — t rn ) is 
the initial Z content at the time of the star formation 
and Q z (m) is the ratio of the mass fraction of new 
metals (or this could refer to some particular clement) 
synthesized in a star of mass m and ejected, to the 
mass of the star. This equation was re-written by 
Maeder (1992) as 

E z (m)=Q z (m) + Z g (t-t m )(R z (m,t)) (21) 

This change is due to a correct use of the definition 
of the yields, which are only the new and not the to- 
tal fraction of metals ejected, so there is no reason to 
subtract Q z (m) a second time. Therefore the chemi- 
cal evolution in our zone can be accounted for in this 
model as 

dZ (t) 1 f* r m d(t'-At) 

^fP = TuJ B(m,t-t')[Q z (m) 

at J g (t) JAtJm d (t') 
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+Z g (t - t')R z (m,t') - R z (m,t')Z g (t)]dmdt' 
+P(t)(Z' g (t) - Z g (t)) (22) 

where Z' g {t) is the metallicity of the neighborhood 
of the zone at time t, and rrid(t') is the mass of the 
star whose lifetime t m is equal to the age of the corre- 
sponding SSP t' at the time at which the computation 
is made. 

Next the code calculates the corresponding metal- 
licity Z g {t) using an equation with a similar form to 
Eq. 19. 

All this is done for each discrete time-step (At) un- 
til the desired time for the region under consideration 
Tq- When t = Tq we obtain the final model results 
for the chemical evolution and star formation history 
of the zone studied. 

The approximation when f g (t) —> 0. To avoid 
any numerical instabilities in the calculation of Z(t) 
when F(t)At ~ f g (t) we use the following approxi- 
mation for Eq. 22: 

^ ft rmd(t' -At) 

Z °® = / B(m,t-t')[Q z (m) 

r (t) J At Jm d (t>) 

+Z g (t - t')R z (m,t')]dmdt' 
+P(t)(Z' g (t) - Z g (t)) (23) 

3.3. The total ejecta 

As explained above, the gas-mass fraction i? z (m, t) 
ejected by a star until a given time t, is obtained from 
the isochrones, but since for a star which does not end 
its evolution as a WD, the core C-ignition phase is the 
last stage included in the above referenced isochrones, 
we do not know the real mass of the remnant. This 
implies that we cannot calculate the total gas-mass 
fraction ejected by a star during the whole of its life 
Rz{m, t m )- We also need the stellar yields, Q z (m), to 
account for the chemical evolution using Eq. 22. Both 
quantities depend on the mass of the star and on its 
initial composition but, unfortunately, there is a lack 
of extensive tables giving these numbers, especially 
for non-solar metallicities. We chose to use the results 
of Renzini & Voli (1981) for the intermediate stars 
(1 - 8M Q ) and those of Maeder (1992) for massive 
stars. In these studies tables are available for solar 
metallicity and also for Z = 0.004 and Z = 0.001. 

Because the code follows the enrichment of the ISM 
in terms of Z, instead of individual elements, we have 
computed Q z (m) for the case of intermediate stars 
by adding all the metal contributors contained in the 
tables of Renzini & Voli for the case of r/ = 1/3 and 
a = 0: the parameter of the rate of mass loss in the 



empirical formulation by Reimers (1975) and the rate 
of the mixing-length to the pressure scale height. For 
all these intermediate stars we have used the remnant 
values only for those stars which do not end their lives 
as a WD because, as was pointed out above, these can 
be obtained directly from the isochrones. For mas- 
sive stars we took both quantities from the tables of 
Maeder (1992), interpolating linearly between the val- 
ues given for the two cases of mass-loss rates for stars 
more massive than 20M Q with solar composition. Fi- 
nally, to obtain both quantities at the compositions 
of the published isochrones, linear interpolations were 
made between Z = 0.004 and Z = 0.02 to obtain the 
values for Z = 0.008 in the case of intermediate stars 
while Z = 0.001 and Z = 0.02 were used to obtain the 
values for Z = 0.004 and Z = 0.008 for massive stars. 
To avoid extrapolation, we used values for Z — 0.004 
for intermediate-mass stars with lower metallicities, 
and values for Z = 0.001 for massive stars of these 
metallicities. For the same reason, we took the solar 
values for higher metallicities for both kind of stars. 

All this information is kept in a file in which we 
interpolate linearly to obtain a more finely divided 
distribution of stars according to their initial mass. 
We must point out that we have taken into account 
only the fractions of new metals ejected Q z (m) when 
the stars die, and not in their previous stages. How- 
ever, this is not the case for R z (m, t) which was com- 
puted from the isochrones. This approximation does 
not introduce any appreciable error, since the frac- 
tion of ejected gas coming from living stars during a 
given period period At is much smaller than that of 
more massive dying stars and therefore the fraction 
of metals that is added is even less important. In any 
case, for these stars their total fraction R z (m,t m ) or 
Q z (m) is always taken into account in the following 
period At. 



21 



Fig. 11. — Comparison of the time variation of the 
metal content and of the fractional mass of gas com- 
puted here with those obtained by AY86 (large sym- 
bols) for different parameter sets (v, ji), where v is 
in units of 10~ 4 Myr~ 1 (the solar neighborhood value 
given by these authors is 1.92 in our units). Mod- 
els with equal fj, have the same point type. To be 
consistent with AY86 our results were obtained us- 
ing a unimodal IMF for which m; = 0.05M Q and 
m u = 60Mq. For v — 1.92 and /i = 1.35 the empir- 
ical age-metallicity relation of Twarog (1980) is well 
reproduced. The wiggles that appear when working 
at high SFR regimes are due to numerical instabilities 
which were diminished by using the approximation 
of Eq. 23. This is caused by the virtual absence of 
available gas. However this effect does not have any 
influence on the observed light of the galaxy, because 
at these high SFR most of the stars we observe at the 
present time were created during the first Gyr. 




5 10 15 

KCyr) 



Fig. 12. — Plot of the contributions of each con- 
stituent SSP to the integrated light in V at 15 Gyr. 
We used the same models as in Fig. 11. Values of the 
parameters are shown beside the corresponding 
curves. 
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3.4. The final distribution of the stars and 
the integrated observables 

To obtain the final distribution of the stars we scan 
all the created SSPs, with ages between and Tq, and 
mark those stars whose lifetimes allow them to be still 
emitting light at the present time Tq. The number 
of stars N(m,t) corresponding to each SSP and which 
survive to the present time, with masses in the range 
(m,m + Am) formed in the time interval (t,t — At) is 
calculated using: 

N(m, t) = M G (t)AmAt (24) 

Once this distribution of stars is known we inte- 
grate the fluxes Fx of all surviving stars of all the 
SSPs in each photometric band 

F x = / N(m,t')F Xi {m,t')dmdt' (25) 

J At Jm l 

where m t > is the mass of a star which dies at t = 
t'. Finally, the code calculates the integrated set of 
colors. 

We proceed as in 2.4 to calculate the integrated 
line indices but now taking into acccount that the 
code yields a stellar distribution that is composed of 
the set of SSPs whose stars are still alive. For each 
surviving star of each created SSP, the code takes the 
log g and the T e ff directly from the isochrones and 
generates the metallicity from the chemical evolution 
in Z which is transformed to [M/H], as in Eq. 10, but 
now Z is the metallicity of the closest isochrone and 
not exactly the one given by the model output. A lu- 
minosity class, as defined in Eq. 9, is assigned to each 
contributing star. Then it calculates the individual 
contribution of a star to the desired line index via the 
corresponding fitting functions and integrates: 

w _ lit Q' W(m,t')N(m,t')F c (m,t')dmdt' 
lit C N(m,t')F c (m,t')dmdt> 

(26) 

where F c (m, t') is the flux of the continuum of the star 
of mass m at time t' obtained by a linear interpolation 
as in Eq. 12. 



3.5.1. The free INP U T parameters 

The code needs the following free input parameters 
which will determine the evolution of the region under 
study: 

• The SFR coefficient v (in units of lO^Myr" 1 ). 

• The minimum gas fraction below which no new 
stars can be created, f 9min ■ In the present study 
we will leave f 3min at 0. 

• The shape of the IMF: the unimodal or the bi- 
modal IMF as were defined in Section 2.2. The 
slope li. The lower and upper mass-cutoff: mi 
and m u . In the present study these are fixed 
at 0.0992 and 72 M Q . We shoulf poiny out 
that higher values for the lower mass-cutoff im- 
ply fewer low-mass stars, and by inference more 
high-mass stars so that the chemical evolution 
proceeds more rapidly and higher mctallicities 
are reached. On the other hand a decrease 
in mi implies that a very important fraction 
of the mass is locked into low-mass stars pre- 
venting ISM enrichment, especially for IMF's 
with /i > 1.35. Generally, using a low IMF 
slope (li < 1.35) higher mctallicities are reached 
(more rapidly if using a high star formation co- 
efficient i/), due to the relative importance of 
massive stars. 

• The age To of the zone and the time-step At. 
Because results do not depend strongly on At 
we mainly use a time-step of 100 Myr~ x (see 
Section 3.5.3 for discussion). 



3.5. Input to the model 

In this section we investigate the influence of the 
most important input parameters; basically we focus 
our attention on i>, fi and the age. These are the 
main parameters determining the star formation and 
the chemical evolution of a desired region of a galaxy, 
and therefore by varying them we can obtain model 
output that can be compared with observations. 
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Table 6: The relative errors in the integrated colors and indices obtained by comparing two time-steps: At = 100 
and 50Myr. We selected values of which together with the age do not allow the metallicity to be greater than 
Z=0.1. The given numbers are the difference in the color or index in magnitude, or this difference as a relative 
fraction (divided by the index at At = 100), when the index is expressed in EW (following the convention of W94). 
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1 
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1 
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2 
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2.446 
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965 
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1.696 


1 
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1.703 


1.301 


1.560 


1.691 


1 


255 


1.443 


1.562 


Fe5709 


0.464 
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0.855 


0.746 


0.926 


0.970 
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0.933 


0.956 


0.852 


0.922 


0.946 





811 


0.859 


0.883 


Fe5782 


0.323 





554 


0.699 


0.595 


0.762 


0.817 





655 


0.768 


0.807 


0.682 


0.760 


0.799 





636 


0.692 


0.728 


NaD 


1.278 


1 


943 


2.607 


2.047 


2.952 


3.453 


2 


375 


3.100 


3.514 


2.543 


3.108 


3.501 


2 


549 


2.998 


3.383 


TiOj 


0.029 





031 


0.033 


0.029 


0.037 


0.041 





031 


0.039 


0.042 


0.034 


0.040 


0.042 





033 


0.036 


0.039 


TiO TI 


0.032 





039 


0.047 


0.037 


0.057 


0.067 





043 


0.061 


0.068 


0.049 


0.062 


0.068 





047 


0.056 


0.062 


Cajjl 


1.151 


1 


604 


1.689 


1.644 


1.718 


1.702 


1 


689 


1.686 


1.679 


1.668 


1.654 


1.657 


1 


622 


1.612 


1.626 


Ca z/ 2 


3.244 


4 


204 


4.339 


4.271 


4.353 


4.256 


4 


337 


4.252 


4.179 


4.274 


4.174 


4.126 


4 


133 


4.041 


4.016 




2.760 


3 


669 


3.772 


3.773 


3.747 


3.580 


3 


812 


3.618 


3.489 


3.723 


3.539 


3.443 


3 


601 


3.433 


3.354 


M ST 


0.449 





643 


0.765 


0.681 


0.820 


0.854 





746 


0.818 


0.844 


0.760 


0.809 


0.837 





745 


0.780 


0.816 



Table 7: The model obscrvablcs obtained with our full chemical evolutionary model for a constant unimodal and 
bimodal IMF with slope 1.35. The SFR coefficient v is in units of 10~ i Myr~ 1 , while the age is in Gyr. Z en d 
indicates the metallicity obtained at the assumed age while < Z > indicates the average metallicity 
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3.5.2. Effects of changing fj, and v 

The main parameters determining the evolution of 
the metallicity in a model are the coefficient v of the 
SFR and the shape of the IMF, mainly its slope (i. 
These two as well as the assumed age for the zone 
under study will determine the final metallicity and 
the average value. 

We have compared our evolution of the metallic- 
ity and the fractional gas mass with AY86, using the 
same IMF (with the same values of the lower and up- 
per mass limits), a linear dependence of the SFR on 
the fractional gas mass, and without using any star 
formation threshold. The evolution of the metallicity 
of the gas, and the gas fraction are plotted in Fig. 11 
for different pairs of model parameters De- 
spite the use of updated isochrones and stellar yields, 
the latter for the case of massive stars, the agree- 
ment is very good. Here also the empirical relation 
between age and metallicity of nearby stars (Twarog 
1980) is fitted for v = 1.92 (the solar neighborhood 
value given in AY86 but in our units of 10 _4 Myr _1 ) 
and a Salpeter IMF (pi — 1.35) . Despite the simi- 
larity of the chemical evolution, their resulting colors 
arc different from ours. This is caused by the use of 
different isochrones and different isochrone-color con- 
versions. This point was already discussed in Bressan 
et al. (1994), who used the same set of isochrones we 
used, but different conversions. 

This figure can be used to see the effect of vary- 
ing the SFR coefficient together with the IMF slope. 
When using a low v only a small fraction of the avail- 
able gas at each time t will be used to form new stars, 
and therefore the ISM enrichment is very slow, espe- 
cially for a high fi, which favors low mass stars which 
do not eject large amounts of matter. In Fig. 12 
where we have represented the contribution of each 
SSP to the presently observed light in the V band, 
we see how most of the observed light comes from 
the recent SSPs. However, for high SFR models the 
amount of available gas rapidly decreases (see the case 
of v = 19.2 in Fig. 11), especially when using high 
IMF slopes, causing the matter to be locked up in 
low-mass stars. In that case the observed light comes 
mostly from the oldest SSPs. For high SFR the metal- 
licity usually increases much faster, but for v > 10, at 
a certain time the metallicity reaches a maximum and 
starts to decline due to ejection of unprocessed mat- 
ter from the oldest and numerous low-metallic SSPs 
stars. Obviously, here one should understand that the 
highest resulting metallicity is reached when combin- 
ing strong i/'s with low /it's. This effect can also be 
seen when combining a normal v = 1.92 with /j, = 0.35 
obtaining Z > 0.1 values (see Fig. 11). At Z = 0.1 



the models are still not good enough to be compared 
with observations. 

Following Arimoto & Yoshii (1987), a typical time- 
scale value for elliptical galaxies is in around ^ 100 x 
10~ 4 Myr~ 1 . However, to fit different kinds of ob- 
jects, often other values over a wide range (0.1 — 200) 
are needed. For the same reason in practice, a con- 
siderable range in \i around 1.35 (Salpeter) might be 
needed in order to fit the observations. 

3.5.3. The time-step 

At is the time-step used in the calculation. In prin- 
ciple this parameter should be infinitely small but one 
has to find a value sufficiently large that the model 
does not spend too much computing time, and suf- 
ficiently small that the results remain accurate. In 
theory, this parameter has to be smaller than the 
cooling time of a cloud to form a new SSP. Empir- 
ically we have made tests to ascertain for which At 
we can still obtain reliable results. Table 6 shows the 
dependence of the final results on this parameter for 
a reasonable range of values (we selected values of 
v and jU which at the assumed age do not allow the 
metallicity to be greater than Z=0.1, see Fig. 11). In 
general the relative error for the colors and indices is 
smaller than a few percent, except for very high SFR 
coefficients, when, for v = 50, this error is around 
<~ 20%. The worst colors again are U-V and V — K . 
This comparison shows that for most cases, a value of 
100 Mgr. is reasonable and is sufficient to obtain rea- 
sonable results. However, when working with high v 
and/or low \x it is advisable to use shorter time-steps. 
Obviously, the use of much shorter At is more impor- 
tant for the study of systems with recent starbursts, 
but is not necessary for old systems. In fact, a useful 
compromise can be to use shorter time-steps for the 
last 100 Mgr. time interval only. 

3.5.4- The static option 

If we remove all the evolutionary aspects of our 
code we obtain a model consisting of a unique SSP as 
in Section 2. This option could in fact be obtained by 
equating the time-step At to the desired age of the 
galaxy, and by using the initial metallicity through- 
out. 

3.6. The observables from our closed-box evo- 
lutionary model 

Here we assumes that our region does not inter- 
change any matter with the neighborhood, the so- 
called closed-box model. In Table 7 we present the 
observables and metallicities (the final and the av- 
erage values) obtained with unimodal and bimodal 
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IMFs constant in time, with slope fi—1.35 for differ- 
ent SFR coefficients and ages. For the sake of brevity 
and since we want to apply our model to early-type 
galaxies we present in this Table results obtained only 
for ^=1.35 because with this value we obtain metallic- 
ities which are not too far from solar, as can be seen in 
Fig. 11. As expected, the bimodal IMF always leads 
to higher metallicitics since the weight of low-mass 
stars is diminished and therefore the chemical evolu- 
tion is more influenced by massive stars. However, as 
explained above, for real galaxies one might need a 
range in [i. This Table and others with further data 
can be obtained from the authors. 

4. Application of our model to three early- 
type galaxies 

In this section we apply the model we described 
in the previous sections to the central regions of a 
few standard early- type galaxies. We first apply the 
static version of our models, and later the evolution- 
ary version. Although static studies have been carried 
out before (e.g. Peletier 1989, Worthcy et al. 1992, 
etc.) in this paper and its companion (Paper II) many 
more spectral line indices are analyzed than is gener- 
ally done. The evolutionary aspect of our model is 
completely new. Up to now, nobody has tried to fit 
at the same time colors and absorption line strengths 
of elliptical galaxies with an evolutionary model. For 
that purpose we have obtained high quality spectro- 
scopic observational data of two standard giant el- 
lipticals (NGC 3379 and NGC 4472) and the bulge 
of the Sombrero Galaxy (NGC 4594). The details 
of these observations (including the conversion of the 
line-strengths to the Lick system) are to be found in 
Paper II. Here as a sample we will analyze the values 
for the three galaxies at 5 arcsec from the center (to 
avoid problems related to the seeing): along the mi- 
nor axis for the Sombrero galaxy and along the major 
axis for the giant ellipticals. The colors for the two 
giant ellipticals are taken from Peletier et al. (1990a) 
and Peletier et al. (1990b) while those for the Som- 
brero galaxy are taken from Hes & Peletier (1993). In 
this paper we will concentrate on what we can learn 
about the models, while in Paper II the observations 
are discussed and emphasis is put on what we can 
learn about the stellar content of elliptical galaxies. 

4.1. Applying the single-age stellar popula- 
tion model 

In Figs. 13 and 14 we have plotted the models, 
together with our three galaxies, in two important di- 
agnostic diagrams: B — V vs. Mg 2 , and of V — K vs. 
Mg 2 . In the first place, one sees that for each galaxy 
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Fig. 13. — Synthetic diagram of B — V vs. Mg2 
obtained by calculating single-age stellar populations 
(SSP) for different ages (ages range from 1 to 17 Gyr. 
along each curve shown) and metallicities, for uni- 
modal IMFs with slopes: 0.35, 1.35 (Salpeter case) 
and 2.35. Also indicated are the observed values for 
the three galaxies at 5 arcsec from their centers. Here 
it is possible to see that a higher than (or at least 
equal to) solar metallicity is required to fit real data, 
whatever the age or the IMF slope. This confirms the 
results given in W94 for a Salpeter IMF. 
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Table 8: Sets of models with parameters (v,l-io) which 
drive the chemical evolution to reach values ~ 2.5 
times solar, and solar metallicity at times, to, of 0.2, 
0.5, 1 and 2 Gyr. 
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there are models which fit the data. In the two dia- 
grams our models seem to fit better than, for example, 
those by W94. Looking at the two diagrams one can 
also see that solar metallicities or larger values are 
required to fit the observational data for the three 
galaxies, whatever their ages or their IMF slopes. 
This result was inferred by ourselves (Casuso et al. 
1996) from an analysis of the Mg2 index alone. 

Comparing models by various authors (see Sec- 
tion 2.6) we have seen that V — K is a very diffi- 
cult color to model, due to uncertainties in the color- 
temperature and color-color transformations for low- 
mass stars and in our knowledge of advanced stages of 
stellar evolution. On the other hand it is also a very 
sensitive color, due to its large wavelength baseline. 
In Fig. 14 we see that V — K depends considerably on 
the IMF slope /i. We consider IMF slopes of 1.35 and 
2.35 and the two shapes of IMF proposed in Section 
2.2. In this plot we see that low-mass stars are impor- 
tant in determining V — K, since there is a substantial 
difference between the predictions for the unimodal 
and bimodal IMFs, especially for high IMFs slopes 
and greater ages. Therefore the stars with masses be- 
low 0.6M Q have a non negligible effect, and can make 
the V — K color redder by ~ 0.4 mag for old popula- 
tions (see for details Section 2.7.1). 

Finally we see that if one includes a small amount 
of extinction, which reddens the colors but does not 
change the line indices, the quality of the fits de- 
grades. In Paper II we will discuss the fit to the other 
colors and absorption lines in more detail. 

4.2. Applying the chemo-evolutionary model 

In an exploratory attempt to find an acceptable 
evolutionary fit to our galaxies we will focus our at- 
tention on various very important observables: the 
colors V — K and U — V, and the Mg 2 and H/3 indices. 
We first run models varying the two main parameters 
(/U,^) and also the ages assigned to the galaxies. Here 
we also study the two forms of the IMF: the unimodal 
and the bimodal cases. In all these cases we use an 
IMF slope constant in time in a closed-box approxi- 
mation. Figs. 15, 16 and 17 show the model output 
for the selected observables when we sweep our pa- 
rameter space. In these plots one can easily see how 
our evolutionary scheme of stellar population synthe- 
sis cannot fit the real data. In particular we see that 
the Mg 2 index of the models is always too low for 
given colors. There is an exception when using ex- 
tremely steep unimodal IMFs for which we can obtain 
Mg 2 > 0.28 (independently of the SFR that we use) 
but which cannot increase as rapidly as the V — K 
color, due to the coolest very low MS stars. 
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Fig. 14. — Shown are various SSP models together 
with observations for our three galaxies. Also shown 
are models from W94 for comparison. 
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Fig. 15. — Synthetic Mg2-(V — K) diagram for four 
SFR regimes (v=l, 5, 20 and lOOx lO^Myr" 1 ), for 
a unimodal IMF with different slopes, from to 4 
in steps of 0.5. The age for the galaxies was varied 
from 1 to 17 Gyr, which corresponds to the last point 
on each curve (where the IMF slope is annotated). 
Observations for the three galaxies are shown. 
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Fig. 16. — The H/3-(U-V) diagram corresponding to 
the same set of parameters as in Fig. 15. 
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Fig. 17. — The (V - K)-Mg 2 diagram obtained for 
the case of a bimodal IMF. 



It is didactically instructive to compare the results 
of our models with the distribution of stars as a func- 
tion of the metallicity with those resulting from ana- 
lytical models. In Fig. 18 we show the contributions 
at the present time of the stars as a function of metal- 
licity and age for some of the best fitting (although 
not acceptable) models, inferred from Figs. 15, 16 and 
17. We have also plotted some of the so-called sim- 
ple analytical models which assume a closed system, 
initially all metal-free gas and a constant IMF (see 
the review of Audouze & Tinsley 1976). We see that 
both, the simple models and ours, yield distributions 
which contain important fractions of low-metallicity 
stars, explaining the mismatch in the V — K vs. Mg2. 
We are clearly dealing with a version of the G-dwarf 
problem in which we have to account for the absence 
of low metallicity stars. This is why, with our invari- 
ant IMF models, we cannot obtain higher values for 
Mg2. If we choose a low value for v to decrease the 
importance of the first SSPs, then we get more or less 
continuous star formation, and the youngest SSPs will 
be the brightest, and since they are young, their Mg2 
will be rather low, even if they are metal-rich. All 
this is thoroughly discussed in Casuso et al. (1996) 
in which we focused our attention exclusively on the 
Mg 2 index. If one looks at the results presented in 
AY86, one sees that the authors have the same prob- 
lem, and that they can barely find a solution that fits 
a typical B — V and V — K for a galaxy. They would 
not have been able to fit its Mg 2 index with their 
model. The previous conclusions are more solid if we 
look at the H/3 vs. U — V diagram, which cannot be 
fitted cither. Finally, in Figs. 15 and 17 there are a 
few models, with very low IMf slopes, which almost fit 
the V — K vs. Mg 2 diagrams of this set of galaxies. 
These models however do not fit the corresponding 
H/3 vs. U — V diagrams (see Fig. 16). 

The fact that the metallicity distribution of our 
model agrees with the predictions of a simple ana- 
lytic model does not mean however that the latter are 
adequate for handling in detail the variations which 
have in fact occurred in stellar populations, even in 
less complex systems such as ellipticals. The present 
model is much more flexible in including the evolu- 
tionary histories of stars (such as mass loss) and al- 
lows us to understand better the star formation his- 
tory. The fact that the stellar frequency distribution 
with metallicity can be reproduced with an analytical 
model does however give some support to our proce- 
dure. 

4-2.1. A variable IMF scenario 

The key question is now how to build a scenario 
that produces a dominant old but metal-rich popu- 



29 



Fig. 18. — Histogram of the cumulative mass fraction 
of living stars as a function of the metallicity for two 
models, with constant IMF, which give fits that are 
the best, but not acceptable, compared to real data. 
'U' means unimodal IMF while 'B' means bimodal. 
These models were selected by looking at Figs. 15, 
16 and 17. We have also plotted a couple of simple 
analytical models fixing the final metallicity to match 
ours, but varying the net yield 7. To be compared 
with these analytical models, in the embedded figure 
(same scale) we plot one of our models that give rise 
to chemical evolution which yields to solar metallic- 
ity. Notice that the analytical models almost match 
ours. All our models assume an age of 17 Gyr. In the 
second figure we plot the corresponding chemical evo- 
lution for these two models. Finally, in the last figure 
we show I re ;(V), the fraction (not cumulative) of the 
total luminosity per time step in V, at the present 
time, of the SSPs formed at different epochs. Notice 
that the cumulative fraction of stars with metallici- 
ties lower than solar is in the range 8 — 36%. Also 
notice that these stars contribute substantially to the 
present light, as indicated in the last plot. 



lation in our observed galaxies. To tackle this prob- 
lem, in this paper we introduce a time dependence 
in the IMF slope as one of the ways which are able 
to yield a metallicity distribution which is skewed to- 
ward higher metallicities. Another would be to in- 
clude some gaseous infall, but we considered this to 
be less lively for ellipticals. In the beginning, for a 
period of time (to) of ~ 1 Gyr, or even less, we fa- 
vor the production of massive stars (using an IMF 
with low slope /ito), while later mainly low-mass stars 
are produced (this IMF slope is termed /z). The first 
to propose this variable IMF scenario was Schmidt 
(1963). It also was analyzed by Arnaud et al. (1992), 
Worthey et al. (1992) and more recently by Elbaz et 
al. (1995) among others. Some authors tackled the 
problem by keeping the IMF slope constant (with the 
Salpeter index), but changing the lower-mass cutoff. 
In our scenario, at the beginning short-lived massive 
stars quickly enrich the ISM but these stars, at the 
same time, do not contribute to the light we are ob- 
serving now. This light comes essentially from stars 
created later from metal-rich gas and a steeper IMF. 

This scenario is plausible physically because: 

• It is likely that the central density and temper- 
ature were initially, high, which favored high- 
mass star formation (see starbursts references) 
in these massive galaxies, so that the IMF slope 
could well have been smaller than 1. 

• Later, conditions tend towards those of the solar 
neighborhood, so that \i becomes larger. 

If we do not want to include a specific treatment 
for the galactic equilibrium between SN rates and 
the galactic mass-dependent binding energy to pre- 
dict when gas densities are able to form stars, we have 
to introduce a phenomenological parameter, which is 
the duration of this short period in which massive 
stars were favored. For this purpose we first plot in 
Fig. 19 the chemical evolution during the first 2 Gyr. 
for flatter IMF slopes (-1 < /z < 1) for different SFR 
regimes (characterized by the v coefficients) . The re- 
sults here do not depend very much on whether we 
choose a unimodal or bimodal IMF (in the latter case 
the metallicities that are reached are slightly higher). 
Because we are limited by the isochrones, which are 
available only for solar and 2.5 times solar metallici- 
ties, we will analyze those chemical evolution param- 
eters for which the two metallicities are reached in 
different periods of time: 0.2, 0.5 1 and 2 Gyr. Ta- 
ble 8 summarizes these results. The main difference 
for two pairs of chemical evolution models which reach 
the same metallicity level at this initial period is that 
for lower /zq the contribution of the earliest SSPs to 
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the present light is less important. 

The (V — K)-Mg2 diagram is a useful diagnostic, 
and helps us to check which solutions are reasonable. 
However, other diagrams such as U — V vs. H/3 serve 
the same purpose (see Figs. 15 and 16 or Figs. 21 and 
22). In Figs. 20, 21 and 22 we plot the diagrams for 
the models with v and fiQ which reach a metallicity of 
2.5 times solar at to following Table 8. Focusing our 
attention on the V — K vs. Mg2 plots we clearly see 
that with a bimodal IMF our data can be fitted bet- 
ter. In these plots we also find that better fits can be 
obtained for to < 1 Gyr. The reason for that is clear: 
if we extend the initial period of time we need lower 
SFR coefficients implying that the most recent stel- 
lar populations become relatively more important and 
therefore the Mg2 decreases. The highest values for 
the Mg 2 index are in fact obtained for t = 0.2 Gyr. 
In Fig. 23 we detail the different contributions to the 
present time of some of the most representative fits. 
We see that the contribution of the stars with metal- 
licities lower than solar, formed during the very short 
early stage of the galaxy evolution, is almost negligi- 
ble (lower than ~ 5%). To illustrate this, we have also 
compared these results with a simple analytic model 
which yields the same final metallicity as our variable 
IMF model. In the same figure we also plotted an 
analytical infall model, in which metal-free gas enters 
at a rate assumed to equal the stellar birthrate (Lar- 
son 1972, Audouze & Tinsley 1976 for a review). We 
see that this model yields the same results as our nu- 
merical variable IMF model. Therefore there are two 
main ways to obtain the required metallicity distribu- 
tions. We have chosen a variable IMF scenario rather 
than infall in our numerical models because this would 
seem to be a more reasonable physical assumption for 
ellipticals and bulges, whereas long-term infall would 
clearly be a reasonable assumption for the evolution 
of the disk population such as that of the solar neigh- 
borhood. 

Finally, since the Mg2 should decrease if we choose 
models (v,^) which reach solar metallicity at to we 
also took to = 0.2 Gyr to make Mg2 as high as pos- 
sible. Fig. 24 shows that for these models we cannot 
obtain the observed high values. The Sombrero is the 
only galaxy which seems to be fitted by this level of 
metallicity for a remaining fi lower than 1.7, but this 
is not a real fit, because with these IMF slopes the 
metallicity continues increasing to reach values higher 
than Z—0.1 and even more, and therefore the result- 
ing observables are not fitted. NGC 3379 seems to 
require higher metallicity than solar but not as high 
as 2.5 times solar, while NGC 4472 seems to be well 
fitted by a value of 2.5 times solar. 

Here we should direct our attention to the recent 



paper of Gibson (1996) in which he criticizes this vari- 
able IMF scenario. He constructs a coupled photo- 
metric and chemical evolution package based on the 
predictions for the bimodal star formation/IMF sce- 
nario of Elbaz ef al. (1995) implementing both the 
isochrones of W94 and those of BBCFN for single 
stellar population models. His models produce V — K 
colors ~ 0.5 — > 1.4 magnitudes too red compared to 
data from ellipticals. We have shown in this section 
that we can fit the data very well. In our opinion the 
reason for the discrepancy found by Gibson (1996) is 
not due to a variable IMF scenario but to the follow- 
ing effects: i) our method of converting the isochrone 
parameters to colors gives a bluer V — K, as can be 
seen for example in Fig. 3 for a single stellar popula- 
tion model of 12 Gyr and metallicity higher than so- 
lar, ii) for his calculations Gibson (1996) used a SSP 
while we have followed the chemical evolution during 
the whole history of the galaxy using our evolutionary 
model. 
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Fig. 19. — The chemical evolution in the first 2 Gyr 
for different values for a bimodal IMF. The hor- 
izontal lines indicate solar and 2.5 times solar metal- 
licities. 
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Fig. 20. — The (V — K)-Mg 2 diagram obtained using 
the assumption that during an initial short period 
of time (0.2, 0.5, 1 and 2 Gyr.) the IMF was less 
steep (with an slope of juo) than during the remain- 
ing time (with an slope of /i). A unimodal IMF was 
used, (f ,/io) are given in brackets while fi is shown at 
the end of each curve, which represents the synthetic 
values obtained for assumed ages ranging along the 
curves from 8 to 17 yr. Notice that the scale is dif- 
ferent from that in Figs. 15, 16 and 17. Following 
Table 8, values of (v,no) were selected to obtain a 
metallicity around 2.5 times solar at the end of the 
initial period of time. On the diagrams we have rep- 
resented only the results obtained for specific models 
(f,/iio) which bracket the observed data. Therefore, 
acceptable fits will be contained for (v^o) values be- 
tween these two limiting models. 
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Fig. 21. — The same as in Fig. 20 but for a bimodal 
IMF. Notice that here the models fit better than in 
the previous unimodal case. 
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Fig. 23. — The same plots as in Fig. 23 but now using 
two of the most representative fits of the variable IMF 
models which fit the data, by inspection of Figs. 20, 
21 and 22. For comparison we represent a simple ana- 
lytical model with a final mctallicity of 2.5 Zq. Notice 
that the cumulative fraction of stars with mctallicities 
lower than solar is now below ~ 5%. Also notice that 
the contribution to the present light of the stars born 
at the early stages of the galaxy evolution (with low 
metallicitics) is much less important than in Fig. 18. 
Finally, we see that another possible scenario capable 
of explaining the observed distribution of stars is by 
means of an infall model in which metal-free gas en- 
ters at a rate assumed to equal the stellar birthrate 
(Larson 1972). 



Fig. 22. — The H/3-(U-V) diagram corresponding to 
the same set of parameters of Fig. 21. 
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5. Conclusions 



We have developed a stellar population model to 
apply to early-type galaxies. It produces optical and 
near-infrared colors, and absorption line strengths in 
lines from 4100 - 8800 A on the Lick system, is ap- 
plicable to systems of intermediate or old age, and 
metallicity larger than 0.1 Z©. The model is chemo- 
evolutionary, i.e. it calculates the properties of a 
stellar system, starting from a primordial gas cloud. 
However, it can easily be used to predict the prop- 
erties of systems with a single age and metallicity. 
The model colors and line strengths are determined 
by integrating stellar observables along theoretical 
isochrones, in this way obtaining Single Stellar Popu- 
lations (SSPs), and then integrating these SSPs over 
time. The model uses isochrones with solar metal 
abundance ratios. As far as possible empirical calibra- 
tions have been used to convert theoretical isochrone 
values to observables for individual stars. Some prop- 
erties of the models are: 

• Extensive comparisons with models from other 
authors show that there is broad overall agree- 
ment in the colors and line strengths. We have 
also made independent estimates of the errors 
in our observables. 

• Our chemical evolutionary model is in good 
agreement with Arimoto & Yoshii (1986). We 
present optical and IR colors which are likely to 
be improved as a result of new calibration rela- 
tions. Also, for the first time, we give integrated 
line strengths for an evolutionary model. 

• We find that for a closed box approximation the 
most metal rich elliptical galaxies cannot be fit- 
ted with a single IMF that is constant in time. 
To solve this problem, we propose a scenario in- 
voking an IMF skewed towards high-mass stars 
during a short, initial period (smaller than 1 
Gyr) , followed by preferential low-mass star for- 
mation in the remaining time. 

• We have briefly tested the model here by fitting 
it to a few key colors and line strength indices, 
using a new data set for 3 standard early-type 
galaxies. In general, satisfactory results are re- 
ported. A much more comprehensive compari- 
son of theory with observations is given in Paper 
II. 
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Fig. 24. — The (V — K)-M.g2 diagram for models 
which reach solar metallicity in an initial period of 
0.2 Gyr. for a bimodal IMF. We can see that the Mg 2 
index is always lower than the observational data. As 
can be deduced from Fig. 21 the Mg2 index is ex- 
pected to decrease with increasing to- The metallic- 
ities in the three galaxies must therefore clearly be 
higher than solar. 



To conclude, we need to make a few statements 
about the applicability of this model. Given the na- 
ture of this type studies, the numbers given in the 
text will soon cease to be the most accurate possible, 
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because better isochrones become available, or better 
calibrations linking one parameter to another. There- 
fore, we will try to update the model as time goes on, 
and interested people can always obtain the most re- 
cent version electronically from the authors. There 
are however a few areas in which we think that fur- 
ther effort by the astronomical community is needed 
to improve models of this kind. These are: 

• We show that absorption line strengths are gen- 
erally more accurate than integrated colors. To 
change this situation better color-color relations 
are needed, especially for very high and very low 
metallicities. 

• Inclusion of more absorption line observations 
will make this kind of models more useful. Espe- 
cially in the near-UV or near-IR very little work 
has been done, except in the region of the Ca II 
triplet. More libraries like the one of Worthcy 
et al. (1994) are urgently needed. 

• It looks as if a substantial part of the uncertain- 
ties are due to incorrect treatment of the later 
stages of stellar evolution, such as the AGB. 
Theoretical work in this area would be very 
valuable. 
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